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Abstract 

The effects of synchrotron radiation on particle motion in storage rings are 
discussed. In the absence of radiation, particle motion is symplectic, and the 
beam emittances are conserved. The inclusion of radiation effects in a classical 
approximation leads to emittance damping: expressions for the damping times 
are derived. Then, it is shown that quantum radiation effects lead to excitation 
of the beam emittances. General expressions for the equilibrium longitudinal 
and horizontal (natural) emittances are derived. The impact of lattice design on 
the natural emittance is discussed, with particular attention to the special cases 
of FODO, achromat, and TME style lattices. Finally, the effects of betatron 
coupling and vertical dispersion (generated by magnet alignment and lattice 
tuning errors) on the vertical emittance are considered. 

1 Introduction 

Beam emittance in a storage ring is an important parameter for characterising machine performance. In 
the case of a light source, for example, the brightness of the synchrotron radiation produced by a beam 
of electrons is directly dependent on the horizontal and vertical emittances of the beam and is one of the 
main figures of merit for users. Second generation light sources had natural emittances of order 100 nm. 
Over the years, significant improvements in lattice designs have been achieved (see Fig. [TJ, motivated 
largely by user requirements; third generation light sources now typically aim for natural emittances of 
just a few nanometres. In the case of colliders for high energy physics, one of the main figures of merif 
is fhe luminosify, which is a measure of fhe rafe of parficle collisions. Fower emittances allow smaller 
beam sizes af fhe inferacfion point, leading to higher particle density in the colliding bunches, and higher 
luminosity for the same total number of particles in the beam. 

There are of course ways of improving the brightness of a light source and the luminosity of a col¬ 
lider without reducing the emittances: in both cases, for example, the beam current could be increased. 
However, beam currents are generally limited by collective effects such as impedance-driven instabili¬ 
ties, Touschek scattering, or (for colliders) beam-beam effects. Designing and operating a storage ring 
for maximum performance involves a good understanding and control of effects that impact the beam 
emittances. 

In this note, we shall consider the emittance of electron (and positron) storage rings: because of 
synchrotron radiation effects, lepton storage rings are able to achieve very small emittances (of order 
1 nm horizontal emittance, and less than 10pm vertical emittance). We shall begin in Sectionj^by re¬ 
viewing some of the key features of beam dynamics in the absence of synchrotron radiation. In particular, 
an important property of the dynamics in such cases is that the particle motion is symplectic: this has the 
consequence that the beam emittances (which characterise the phase space volume occupied by the par¬ 
ticles in a beam) are conserved as the beam moves around the storage ring. We shall then show that, in a 
classical approximation, radiation effects lead to damping of the emittances. We shall derive expressions 
for the exponential damping times. Then, we shall discuss how quantum effects of synchrotron radiation 
lead to excitation of the beam emittances. As a result, the emittances of beams in electron (or positron) 
storage rings reach equilibrium values determined by the beam energy and lattice design. 

In Sectionj^we shall apply the expression for the natural emittance derived in Sectionj^to partic¬ 
ular styles of lattice design. In particular, we shall consider FODO, double bend achromat, theoretical 
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Fig. 1 : Natural emittance of a number of synchrotron light sources. Reductions in natural emittance have been 
driven by the need to produce higher radiation beam brightness for users. 


minimum emittance, and multi-bend achromat lattices. Double bend achromats are of particular interest 
for light sources, because they achieve low natural emittance (leading to high brightness) while pro¬ 
viding long, dispersion-free (or low dispersion) straight sections that are ideal locations for insertion 
devices such as undulators or wigglers |[T|. Insertion devices are useful for providing intense beams of 
synchrotron radiation with specific properties. 

In a planar storage ring, the vertical emittance is dominated by alignment and tuning errors, rather 
than by the design of the lattice. In Section we shall discuss how the vertical emittance is related 
to a range of errors, including steering errors, tilt errors on quadrupoles and vertical alignment errors 
on sextupoles. Betatron coupling and vertical dispersion are important features of the dynamics in this 
context, and both will be discussed. Optimisation of a lattice design for a low-emittance storage ring 
will generally involve simulations to characterise the sensitivity of the vertical emittance to different 
types of machine error. For this, techniques are needed for accurate computation of the equilibrium 
emittances from models in which different errors can be included. We shall consider three techniques 
that are widely used for emittance computation, discussing the envelope method in particular in some 
detail. Finally, we shall mention briefly some of the issues associated with operational tuning of a storage 
ring for low-emittance operation. 

2 Beam dynamics with synchrotron radiation 

In this section, we shall review the relevant aspects of beam dynamics needed for understanding the 
effects of synchrotron radiation. Our focus will be on electron (or positron) synchrotron storage rings. 
Initially, we shall neglect radiation effects; then, we shall include the emission of synchrotron radiation 
as a perturbation to the motion of individual particles. This approach is valid if radiation effects are 
relatively weak, which means that the energy lost by a particle through radiation in one synchrotron 
period should be small compared to the total energy of a particle. This is almost invariably the case 
for practical storage rings. We shall consider only incoherent synchrotron radiation, in other words 
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we shall assume that the motion of each particle and the radiation that it produces can be considered 
independently of all other particles in the beam. In some regimes (including, for example, in free electron 
lasers) particles generate radiation coherently, leading to a strong enhancement of the radiation produced 
by a beam. Generally, some special efforts are needed to achieve the generation of coherent synchrotron 
radiation with sufficient intensity that it has a measurable impact on the beam; we shall not discuss such 
situations here. 

Briefly, we shall proceed as follows. The symplectic motion of particles in an accelerator (i.e. mo¬ 
tion neglecting synchrotron radiation and collective effects) is conveniently described using action-angle 
variables. We shall define these variables, and use them to review the key features of particle motion in 
synchrotron storage rings. We shall then include the effects of synchrotron radiation, initially in a classi¬ 
cal approximation, leading to expressions for the energy lost per turn in a storage ring, and the damping 
times for the horizontal, vertical and longitudinal emittances. Finally, we shall discuss the effects of 
quantum excitation, and derive results for the equilibrium beam emittances. These results will be used 
in Section where we consider how the equilibrium emittances are affected by the lattice design in a 
storage ring. 

2.1 Symplectic motion 

We work in a co-ordinate system based on a reference trajectory that we define for our own convenience 
(see Fig. |^. The distance along the reference trajectory is specified by the independent variable s. 
For simplicity, in a planar storage ring, the reference trajectory is generally chosen to be a straight line 
(passing through the centres of all quadrupole and higher-order multipole magnets) everywhere except 
in the dipoles. In the dipoles, the reference trajectory follows the arc of a circle with radius p, such that: 

Bp=—, (1) 

q 

where B is the dipole field, Pq is the reference momentum (i.e. the momentum of particles for which the 
storage ring is designed) and q is the particle charge. Bp is the beam rigidity. 

At any point along the reference trajectory, the position of a particle is specified by the x and y 
co-ordinates in a plane perpendicular to the reference trajectory. We follow the convention in which x is 
the horizontal (transverse) co-ordinate, and y is the vertical co-ordinate. 



Fig. 2: Co-ordinate system in an accelerator beam line. The reference trajectory can be defined arbitrarly, but is 
generally chosen so that it describes the trajectory of a particle with momentum equal to the reference momentum 
Pq. The distance along the reference trajectory is parameterised by the independent variable s. At any point 
along the reference trajectory, the transverse position of a particle is specified by Cartesian co-ordinates in a plane 
perpendicular to the reference trajectory. 
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Fig. 3: Longitudinal co-ordinates in an accelerator beam line. The longitudinal co-ordinate z indicates the time 
that a particle crosses a plane perpendicular to the reference trajectory at a position s along the reference trajectory. 


To describe the motion of a particle, we need to give the components of the momentum of a 
particle, as well as its co-ordinates. In the transverse directions (i.e. in a plane perpendicular to the 
reference trajectory) we use the canonical momenta Q scaled by the reference momentum Pq- 
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Px = 

^0 
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Py = 

^0 


dx . 
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I /t 


( 2 ) 

(3) 


Here, m and q are the mass and charge of the particle, 7 is the relativistic factor for the particle, 
and Ax and Ay are the x and y components respectively of the electromagnetic vector potential. The 
transverse dynamics are described by giving the transverse co-ordinates and momenta as functions of s 
(the distance along the reference trajectory). 

To describe the longitudinal dynamics of a particle, we use a longitudinal co-ordinate z defined 
by: 

z =/3oc(fo - f), (4) 

where /3o is the normalised velocity of a particle with the reference momentum Pq, Iq is the time at 
which the reference particle is at a location s, and t is the time at which the particle of interest arrives 
at this location. Physically, the value of 2 for a particle is approximately equal to the distance along 
the reference trajectory between the given particle and a reference particle travelling along the reference 
trajectory with momentum Pq (see Fig. |^. A positive value for z means that the given particle arrives 
at a particular location at an earlier time than the reference particle, i.e. the given particle is ahead of the 
reference particle. 

The final dynamical variable needed fo describe fhe mofion of a particle is fhe energy of fhe 
particle. Rafher fhan use fhe absolufe energy or momenfum, we use fhe energy deviation S, which 
provides a measure of fhe difference befween fhe energy P of a particle and fhe energy of a parficle wifh 
fhe reference momenfum Pq: 

c 1 1 7 ^ 

Here, 70 is fhe relafivisfic facfor for a parficle wifh momenfum equal fo fhe reference momenfum. A 
parficle wifh momenfum equal fo fhe reference momenfum has <5 = 0 . 

Using fhe above definitions fhe co-ordinafes and momenfa form canonical conjugate pairs: 

{x,Px), iy,Py), {Z,6). ( 6 ) 
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This means that (continuing to neglect radiation and collective effects) the equations of motion for parti¬ 
cles in an accelerator beam line are given by Hamilton’s equations Q, with an appropriate Hamiltonian 
that describes the electromagnetic fields along the beam line. In a linear approximation, the change in the 
values of the variables when a particle moves along a beam line can be represented by a transfer matrix, 


R: 


X \ 


/ X \ 
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= ii’(si;so) • 
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Py 


Py 
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5 / 

S=Sl 

V ^ 4 


(V) 


It is a general property of Hamilton’s equations that the transfer matrix R is symplectic. Mathematically, 
this means that R satisfies the relation: 

R^SR = S, (8) 

where S is the antisymmetric matrix: 


/ 0 1 0 0 0 0 \ 

-1 0 0 0 0 0 

0 0 0 1 0 0 

0 0 -1 0 0 0 

0 0 0 0 0 1 

Vo 0 0 0-10/ 


(9) 


The symplectic condition ([^ imposes strong constraints on the dynamics. Physically, symplectic 
matrices preserve volumes in phase space (this result is sometimes expressed as Liouville’s theorem 0). 
For example, for a linear transformation in one degree of freedom, a particular ellipse in x-px phase 
space will be transformed to an ellipse with (in general) a different shape; but the area of the ellipse will 
remain the same. The number of invariants associated with a linear symplectic transformation is at least 
equal to the number of degrees of freedom in the system. Thus, for motion in three degrees of freedom, 
there are at least three invariants. For particles in a beam in an accelerator beam line, the invariants are 
associated with the emittances. If there is no coupling between the degrees of freedom (so that motion 
in any direction x,y or z is independent of the motion in the other two directions) then we can associate 
an emittance with each of the three co-ordinates, i.e. there is a horizontal emittance, a vertical emittance 
and a longitudinal emittance. We shall give a more formal definition of the emittances shortly. 

Consider a particle moving through a periodic beam line, without coupling (i.e. a beam line with 
no skew quadrupoles or solenoids). After each periodic cell, we can plot the horizontal co-ordinate x 
and momentum px as a point in the horizontal phase space. After passing through many cells, observing 
the particle always at the corresponding locations in successive cells, and assuming that the motion of 
the particle is stable, we find that the points trace out an ellipse in phase space. The shape of the ellipse 
defines the Courant-Snyder parameters |j^ in the beam line at the observation point: see Fig. The 
area of the ellipse is a measure of the amplitude of the oscillations. We define the horizontal action Jx 
of the particle such that the area of the ellipse is equal to ttJx- 

Applying simple geometry to the phase space ellipse, we find that the action (for uncoupled mo¬ 
tion) is related to the Cartesian variables for the particle by: 

‘^Jx = IxX^ + 2axxpx + PxpI, (10) 

where the Courant-Snyder parameters satisfy the relation: 

Mx-al = l. ( 11 ) 
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Fig. 4: Ellipse in phase space defined by plotting the co-ordinate x and conjugate momentum px of a particle after 
each pass through a unit cell in a periodic beam line. The shape of the ellipse is described by the Courant-Snyder 
parameters ax, Px and ^x', the area of the ellipse is ttJx, where Jx is the action variable of the particle. The shape 
of the ellipse changes depending on the chosen starting position within a unit cell; the action remains the same for 
any given particle. 


We define the horizontal angle variable as follows: 


tan (px 


Q PX 

Px Oix ■ 

X 


( 12 ) 


For a particle with a particular action (i.e. on an ellipse with a given area) the angle variable specifies 
the position of the particle around the ellipse. The action-angle variables Q provide an alternative to 
the Cartesian variables for describing the dynamics. Although we have not shown that this is the case, 
the action-angle variables form a canonical conjugate pair: that is, the equations of motion expressed 
in terms of the action-angle variables can be derived from Hamilton’s equations, using an appropriate 
Hamiltonian (determined as before by the electromagnetic fields along the beam line). The advantage 
of using action-angle variables to describe particle motion in an accelerator is that, under symplectic 
transport (i.e. neglecting radiation and collective effects), the action of a particle is constant. We can of 
course define vertical and (in a synchrotron storage ring) longitudinal action-angle variables in the same 
way as we defined the horizontal action-angle variables. 

The expressions for the action ( [T0| ) and the angle ( [T2] ) can be inverted, to give expressions for the 
Cartesian co-ordinate and momentum in terms of Jx and (fix'. 


X 


Px 


\/‘^i^xJx cos (j^x 



(sin (px + ax cos (px) ■ 


(13) 

(14) 


The emittance Ex of a bunch of particles can be defined as the average action of all particles in the 

bunch: 

ex = {Jx). (15) 
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For uncoupled motion, and assuming that the angle variables of different particles are uncorrelated, it 
follows from (131 and ( [T4| ) that the second order moments of the particle distribution are related to the 
Courant-Snyder parameters and the emittance: 


{xpx) 

{pD 


) 

Oix^X 
'Jx^x ■ 


(16) 

(17) 

(18) 


Using we then find that the emittance can be expressed in terms of the second order moments as: 


= {x^){pI) - {xpxf- 


(19) 


However, we stress that this relation holds only for uncoupled motion. The expression for the emittance 
( [T5] ) can be generalised without too much difficulty to coupled motion (see, for example 0), leading to 
normal mode emittances that are conserved under symplectic transport even where coupling is present. 
However, the expression for the emittance ( [T9] ) is less easily generalised to include coupling, and an 
emittance that is defined by ([T9|) will, in general, nol be consfanf in a beam line where fhere is coupling. 


2.2 Vertical damping by synchrotron radiation 

So far, we have considered only symplectic transport, i.e. motion of a particle in drift spaces or in the 
electromagnetic fields of dipoles, quadrupoles, rf cavities etc. without any radiation. However, we know 
that a charged particle moving through an electromagnetic field will (in general) undergo acceleration, 
and a charged particle undergoing acceleration will radiate energy in the form of electromagnetic waves. 
We now address the question of the impact that this radiation will have on the motion of a particle in a 
synchrotron storage ring. We shall consider first the case of uncoupled vertical motion: for a particle in 
a storage ring, this turns out to be the simplest case. Since we are primarily interested in the dynamics 
of the particles generating the radiation, we quote a number of results regarding the properties of the 
radiation itself (rather than derive these results from first principles). 

The first result that we quote for the properties of synchrotron radiation, is that radiation from a 
relativistic charged particle is emitted within a cone of opening angle of I/ 7 , where 7 is the relativistic 
factor for the particle Q. The axis of the cone is tangent to the trajectory of the particle at the point 
where the radiation is emitted. For an ultra-relativistic particle, 7 S> 1, and we can assume that the 
radiation is emitted directly along the instantaneous direction of motion of the particle. 

Consider a particle with initial momentum P « Po, that emits radiation carrying momentum dP. 
The momentum of the particle after emitting radiation is: 

P' = P - dP ^ P (^1- . (20) 

Since there is no change in direction of the particle, the vertical component of the momentum must scale 
in the same way as the total momentum of the particle: 

Py 

Now we substitute this into the expression for the vertical betatron action (valid for uncoupled motion): 

2jy = ’^yH + “^OtyPPy + PyPyj (22) 

to find the change in the action resulting from the emission of radiation: 

dP 

dJy = - {ayVPy + jdypl) — . (23) 
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Note that in (231 we neglect a term that is second order in dP/P q. This term vanishes in the classical 
approximation when we consider the emission of an infinitesimal amount of radiation in an infinitesimal 
time interval dt; however, we shall see later that including quantum effects, the second order term will 
lead to excitation of the action. Retaining for the present only the first order term in dP/Po, averaging 
(231 over all particles in the beam gives: 


(^dJy'^ - dSy - Sy 


dP_ 


(24) 


where we have used: 


{ypy) = -aySy, 

(pI) = 

and: 

f^y'yy ~ ^y ~ 


(25) 

(26) 

(27) 


The emittance is conserved under symplectic transport, so if the effects of radiation are ‘slow’ 
(i.e. the rate of change of energy from radiation is small compared to the total energy of a particle 
divided by the revolution period), then for a particle in a storage ring we can average the momentum loss 
around the ring. From ([24]): 


d£y _ ^y f dP ^ Uq _ 2 

dt Tq J Pq EqTq ^ Ty 


(28) 


where Tq is the revolution period, and Uq is the energy lost through synchrotron radiation in one turn. 
The approximation is valid for an ultra-relativistic particle, which has Eq Pqc. The damping time Ty 
is defined by: 

Ty = 2^Tq. (29) 

The evolution of the emittance is given by: 


^yit) 


£y{t = 0 ) exp 



(30) 


Typically, in an electron storage ring, the damping time is of order several tens of milliseconds, while 
the revolution period is of the order of a microsecond. In such a case, radiation effects are indeed slow 
compared to the revolution frequency. 

Note that we made the assumption that the momentum of the particle was close to the reference 
momentum, i.e. P « Pq. If the particle continues to radiate without any restoration of energy, we will 
reach a point where this assumption is no longer valid. However, electron storage rings contain rf cavities 
to restore the energy lost through synchrotron radiation. For a thorough analysis of synchrotron radiation 
effects on the vertical motion (at least, with a classical model for the radiation), we should consider the 
change in momentum of a particle as it moves through an rf cavity. However, in general, rf cavities 
are designed to provide a longitudinal electric field. This means that particles experience a change in 
longitudinal momentum as they pass through a cavity, without any change in transverse momentum. In 
other words, the vertical momentum py of a particle will remain constant as the particle moves through 
an rf cavity, which will therefore have no effect on the emittance of the beam. 

To complete our calculation of the vertical damping time, we need to find the energy lost by a 
particle through synchrotron radiation on each turn through the storage ring. At this point, we quote a 
second result from the theory of synchrotron radiation: the radiation power from a relativistic particle 
following a circular trajectory of radius p is given by Lienard’s formula |j^: 


P,= 


q^c /3V 
Gvreo 


C^c c^P^ 

27r 


2-k ^ 271 ^ 


(31) 







where the particle has charge q, velocity /3c ps c, energy E = 'ymc^ and momentum P = j5^mc. The 
particle travels on a path with radius p in a magnetic field of strength B. The approximation in the final 
expression of ( |3T] ) is valid for ultra-relativistic particles, 7 1 . cq is the permittivity of free space, and 

is a physical constant given by: 


E = ^ 

3eo(mc2)^’ 


(32) 


For electrons, k, 8.846 x 10“^ m/GeV'^. Note that the radiation power has a very strong scaling 
with the particle mass: the larger the mass of the particle, the smaller the amount of radiation emitted. 
In proton storage rings, except at extremely high energy, synchrotron radiation effects are generally 
negligible. For a particle with the reference energy, travelling close to the speed of light along the 
reference trajectory, we can find the energy loss by integrating the radiation power around the ring: 


Uo = 

Using the expression ( [3T] ) for P^, we find: 

Uo 


Pry dt 




p ^ 
^7 ^ • 


^ ds, 


(33) 


(34) 


where p is the radius of curvature of the particle trajectory, and we assume that the particle energy is 
equal to the reference energy Eq. For convenience, we assume that the closed orbit is the same as the 
reference trajectory for a particle with the reference momentum. 

Following convention, we define the second synchrotron radiation integral, I 2 @ : 

h = ^ ds. (35) 

In the ultra-relativistic limit, the energy loss per turn Uq is written in terms of I 2 as: 

Uo = ^Ep2. (36) 

ZTT 

Note that I 2 is a property of the lattice (actually, a property of the reference trajectory), and does not 
depend on the properties of the beam. Conventionally, there are five synchrotron radiation integrals 
used to express the effects of synchrotron radiation on the dynamics of ultra-relativistic particles in an 
accelerator. The first synchrotron radiation integral is not, however, directly related to the radiation 
effects. It is defined as: 

I^ = ds, (37) 

where px is the horizontal dispersion. Ji is related to the momentum compaction factor ap, which plays 
an important role in the longitudinal dynamics, and describes the change in the length of the closed orbit 
with respect to particle energy: 

AC . . 


—— ap6 0{6 ). 

Go 


(38) 


The length of the closed orbit changes with energy because of dispersion in regions where the reference 
trajectory has some curvature (see Fig.|^: 


dC = {p + x) d9 = 


1 + - 
P 


ds. 


If X = pxd, then: 


dC = 


P 


The momentum compaction factor can be written: 

1 dC 


Co d5 


5=0 


Co 


ds. 


Vx 


h 


= ^ f —ds= 


Co 


(39) 


(40) 


(41) 
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Fig. 5; Change in path length of a particle following a trajectory offset from the reference trajectory. If a particle 
has co-ordinate x and follows a path parallel to the reference trajectory, then the length of the path followed by 
the particle is dC = (1 -f x/p)ds, where p is the radius of curvature of the reference trajectory, and ds is the 
corresponding distance along the reference trajectory. 


2.3 Horizontal damping 

Analysis of the effect of synchrotron radiation on the vertical emittance was relatively straightforward. 
When we consider the horizontal emittance, there are three complications that we need to address. First, 
the horizontal motion of a particle is often strongly coupled to the longitudinal motion. We cannot treat 
the horizontal motion without also considering (to some extent) the longitudinal motion. Second, where 
the reference trajectory is curved (usually, in dipoles), the length of the path taken by a particle depends 
on the horizontal co-ordinate with respect to the reference trajectory. This can be a significant effect 
since dipoles inevitably generate dispersion (a variation of the orbit with respect to changes in particle 
energy), so the length of the path taken by a particle through a dipole will depend on its energy. Finally, 
dipole magnets are sometimes built with a gradient, in which case the vertical field seen by a particle in 
a dipole will depend on fhe horizonfal co-ordinafe of fhe particle. 

Coupling befween fransverse and longifudinal planes in a beam line is usually represenfed by fhe 
dispersion, rjx and rjpx, defined by: 


rjx 


Vpx 


dxc 


dS 


5=0 


dp X,CO 

d6 


<5=0 


(42) 

(43) 


where Xco and Px,co are fhe co-ordinafe and momenfum for a parficle wifh energy deviafion 5 on a 
closed orbif. We use fhe horizonfal action-angle variables Jx and (f)x to describe the horizontal betatron 
oscillations of a particle with respect to the dispersive closed orbit, i.e. the closed orbit for a particle with 
energy deviation 5. In terms of the horizontal dispersion and betatron action, the horizontal co-ordinate 
and momentum of a particle are given by: 


X 


Px 


^2PxJx cos (jyx 4“ ’dxd'! 



(sin (f)x + ax cos (l)x) + 'q^xd- 


(44) 

(45) 


When a particle emits radiation, we have to take into account both the change in momentum of the 
particle, and the change in co-ordinate and momentum with respect to the new (dispersive) closed orbit. 
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Note that when we analysed the vertieal motion, we assumed that there was no vertieal dispersion. This 
is the ease in an ideal, planar storage ring, but as we shall diseuss later, alignment errors on the magnets 
ean lead to the generation of some vertieal dispersion that depends on the errors, the effeets of whieh 
eannot always be negleeted. 

Taking all the above effeets into aeeount for the horizontal motion, we ean proeeed along the 
same lines as for the analysis of the vertieal emittanee. That is, we first write down the ehanges in eo- 
ordinate x and momentum resulting from an emission of radiation with momentum dP (taking into 
aeeount the additional effeets of dispersion). Then, we substitute expressions for the new eo-ordinate and 
momentum into the expression for the horizontal betatron aetion, to find fhe ehange in aefion resulting 
from fhe radiation emission. Averaging over all parfieles in fhe beam gives fhe ehange in fhe emittanee 
fhaf resulfs from radiation emission from eaeh parfiele in fhe beam. Finally, we infegrafe around fhe ring 
(faking aeeounf of ehanges in pafh lengfh and field sfrengfh wifh fhe horizonfal posifion in fhe bends) fo 
find fhe ehange in emittanee over one furn. 

Filling in fhe sfeps in fhis ealeulafion, we proeeed as follows. Firsf, we nofe fhaf, in fhe presenee 
of dispersion, fhe aetion is written: 

2Ja; = 7x5^ + 2axXpx + /3 xpI, (46) 

where x and px are fhe horizonfal eo-ordinafe and momenfum wifh respeef fo fhe dispersive elosed orbif: 


X = X - TJxd, 

Px — Px 'dpx^- 

Afler emission of radiafion earrying momenfum dP, fhe variables ehange by: 


<5 

1 -^ 

5- 

Po 

G 

X 

I-)- 

X + rjx- 
1 

Px 

!-)■ 

Px ( 

1 - 


dP 


dP\ dP 

We wrife fhe resulfing ehange in fhe aetion as: 

Jx * y Jx T dJx- 

The ehange in fhe horizonfal aefion is: 

dP /dPV 


dJx = -wi— + ru2 ^ 


P 


Po 


where, in fhe limif —)■ 0: 


mi = axxpx + PxpI - r]x{'yxX + axPx) - VpxiotxX + fixPx), 


(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 


and: 


^2 - „ ipx'Hx P “^(^xPxVpX T f^xPpx) ipxdx T Pxdpx) Px T 


(55) 


Treafing radiafion as a elassieal phenomenon, we ean fake fhe limif dP —?■ 0 in fhe limif of small lime 
inlerval, dt —)■ 0. In fhis approximation, fhe ferm fhaf is seeond order in dP vanishes, and we ean wrife 
for fhe rafe of ehange of fhe aefion: 


dJx 

dt 


1 dP 


7 


(56) 
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where is the rate of energy loss of the particle through synchrotron radiation (jSTj). To find the average 

rate of change of horizontal action, we integrate over one revolution period: 


dJx 

dt 


1 

% 


P'1 , 

wi-^ dt. 
Poc 


(57) 


It is more convenient, given a particular lattice design, to integrate over the circumference of the ring, 
rather than over one revolution period. However, we have to be careful changing the variable of integra¬ 
tion (from time t to distance s) where the reference trajectory is curved: 


dt = 


c 



ds 

c 


(58) 


So: 


dJx 

dt 


1 

ToPoc^ 



(59) 


where the rate of energy loss P^ is given by 

We have to take into account the fact that in general, the field sfrengfh in a dipole can vary wifh 
posifion. To firsl order in x we can write: 


dB 

B = Bo + X—. 

ox 


(60) 


Subsfilufing (601 info ( [3T] ), and wifh fhe use of (54i, we find (after some algebra) fhaf, averaging over all 
particles in fhe beam: 


X 


^ 1^7 ( 1 + - 


ds = cC/o ( 1 - ^ ) Ex, 


(61) 


where fhe energy loss per furn Uq is given by ( [36] ), (he second synchrofron radiation infegral I 2 is given 

2 + 2 A:i I ds. (62) 


by (351, and fhe fourfh synchrofron radiation infegral is 1^: 

1 


h = 


rh 

P 


ki is (he normalised quadrupole gradienf in fhe dipole field: 


ki 


q dBy 

Pq dx 


(63) 


Nofe fhaf in (62 1 , fhe dispersion and quadrupole gradienf confribufe fo fhe infegral only in (he dipoles: in 
ofher parfs of (he ring, where fhe beam follows a sfraighf pafh, (he curvafure l/p is zero. 

Averaging (59) over all particles in (he beam and combining wifh (|6T]) we have: 


d£x _ __ h\ 

dt ~ ToEo\ 


(64) 


Defining fhe horizonfal damping fime Tx'. 


where: 


fx — . TT -fg, 

Jx Uo 


h 

jx — ^ 


fhe evolufion of fhe horizonfal emiffance can be written: 

d£r 2 


dt Tx ^ 


(65) 

( 66 ) 


(67) 
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The quantity jx is called the horizontal damping partition number. For most synchrotron storage ring 
lattices, if there is no gradient in the dipoles then jx is very close to 1. From ( [67] ) the horizontal emittance 
decays exponentially: 

ex(t) = ea;(f = 0)exp^-2^^ . (68) 


2.4 Longitudinal damping 

So far we have considered only the effects of synchrotron radiation on the transverse motion, but there 
are also effects on the longitudinal motion. Generally, synchrotron oscillations are treated differently 
from betatron oscillations because in one revolution of a typical storage ring, particles complete many 
betatron oscillations but only a fraction of a synchrotron oscillation. In other words, the betatron tunes 
are vp ^ 1, but the synchrotron tune is Vg <C 1. To find the effects of radiation on synchrotron motion, 
we proceed as follows. We first write down the equations of motion (for the dynamical variables z and 5) 
for a particle performing synchrotron motion, including the radiation energy loss. Then, we express the 
energy loss per turn as a function of the energy deviation of the particle. This introduces a damping term 
into the equations of motion. Finally, solving the equations of motion gives synchrotron oscillations (as 
expected) with amplitude that decays exponentially. 

The changes in energy deviation 6 and longitudinal co-ordinate 2 ; for a particle in one turn around 
a storage ring are given by: 


Az 


OpCqS, 


cjrfz \ U 


(69) 

(70) 


where V^f is the rf voltage, the rf frequency, Eq is the reference energy of the beam, (/)s is the nominal 
rf phase, and U (which may be different from Uq) is the energy lost by the particle through synchrotron 
radiation. Strictly speaking, since the longitudinal co-ordinate 2 ; is a measure of the time at which a 
particle arrives at a particular location in the ring, changes in z with respect to energy should be written 
in terms of the phase slip factor rjp, which describes the change in revolution period with respect to 
changes in energy, rather than in terms of the momentum compaction factor ap. The phase slip factor 
and the momentum compaction factor are related by (see, for example ||7|): 

r]p = ap-\, (71) 

7 o 

where 70 is the relativistic factor for a particle with the reference momentum. But for a storage ring 
operating a long way above transition (which is the situation we shall assume here) ap » l/ 7 g, so 
Vp ~ Oip- It is slightly more convenient to work with the momentum compaction factor, since this 
depends (essentially) on just the geometry of the lattice and the optical functions (in particular, the 
dispersion); whereas the phase slip factor depends also on the beam energy. 

If the revolution period in the storage ring is Tq, then we can write the longitudinal equations of 
motion for the particle: 


(M 

dt 

dz 

dt 


c ) EoTo' 

—apc6. 


(72) 

(73) 


To solve these equations, we have to make some assumptions. First, we assume that 2 ; is small compared 
to the rf wavelength: 


UJrf\z 

C 


< 1 . 


(74) 
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The synchrotron radiation power produced by a particle depends on the energy of the particle. We assume 
that the energy deviation is small, |(5| <C 1, so we can work to first order in 6: 


U = Uo + AE 


dU 

dE 


E=Eo 


Uo + Eo5 


dU 

Zb 


E=Eo 


(75) 


Finally, we assume that the rf phase (j)s is set so that for z = J = 0, the rf cavity restores exactly 
the amount of energy lost by synchrotron radiation. With these assumptions, the equations of motion 
become: 


ds qVrf Wrf . dU 

dt EqTq c Tq dE 

dz _ 

dt ^ 


E=Eo 


Taking the derivative of (761 with respect to t, and substituting for dz/dt from ([TT]) gives: 


(76) 

(77) 


dH 

dt^ 


+ + iZgS = 0 . 


(78) 


This is the equation for a damped harmonic oscillator, with frequency czg and damping constant a e given 
by: 


(jj 


2 

s 


OiE 


qV^ 

Eo 


(1 

COS((ps) O-p, 

7o 


1 dU 
2To dE 


E=Eo 


If <C iZs, the energy deviation and longitudinal co-ordinate damp as: 

5{t) = <5o exp(-Q;£;f) sm{uJst - 9o), 


(79) 

(80) 


(81) 


ck c 

z{t) = —^ Jo exp(—agf) cos(wsf — 0o)- (82) 

COs 

where Jq is a constant (the amplitude of the oscillation in J at f = 0), and Oq is a fixed phase (the phase 
of the oscillation at f = 0). 

To find an explicit expression for the damping constant oe, we need to know how the energy loss 
per turn U depends on the energy deviation 6. The total energy lost per turn by a particle is found by 
integrating the synchrotron radiation power over one revolution period: 


U = 


dt. 


(83) 


To convert this to an integral over the circumference, we should recall that the path length depends on 
the energy deviation; so a particle with a higher energy takes longer to travel around the lattice: 


dC 


1 


c c 

Therefore, the radiation energy loss per turn is: 

uP- 


dt = — = - lH— ds = - l-|- ds. 


ijxd 


1 + 


Vxd 


ds. 


(84) 


(85) 
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Using (|3T]), we find after some algebra: 


dU 

dE 


= 3z 


■ ^ 


where Uq is given by 


\E=Eo ^0 

and the longitudinal damping partition number is: 

h 

jz = 2 + -^. 


( 86 ) 


(87) 


I 2 and I 4 are the same synchrotron radiation integrals that we saw before, in (351 and (62l. Finally, we 
can write the longitudinal damping time: 


1 2 

- — - 

OtE 3z Uq 


( 88 ) 


Neglecting coupling, the longitudinal emittance can be given by a similar expression to the hori¬ 
zontal and vertical emittance: 

Ez = V {z^){S'‘^) - {zS)'^. (89) 

Even where dispersion is present, so that the horizontal and longitudinal motion are coupled, the expres¬ 
sion ( [89] ) can provide a useful definition of the longitudinal emittance, since the longitudinal variables 
usually have a much weaker dependence on the transverse variables, than the transverse variables have 
on the longitudinal. Since the amplitudes of the synchrotron oscillations decay with time constant Tz, the 
damping of the longitudinal emittance can be written: 

- v) • (90) 


£z{t) = £z{t = 0)exp 


It is worth commenting on the fact that the horizontal, vertical and longitudinal emittances are 
all damped by synchrotron radiation with exponential damping times that depend on the beam energy 
and the rate at which particles lose energy through synchrotron radiation. In the case of the horizontal 
and longitudinal emittances, there is an additional factor in the expressions for the damping times that 
depends on details of the lattice, or, more precisely, on the properties of the dipoles. The additional 
factors are given by the damping partition numbers jx and jz- From ( 66 ) and (|87|), we see that: 


3x + 3z — 3- 


(91) 


In general, there can also be a vertical damping partition number jy, although in the simple case we have 
considered here (of a perfectly planar storage ring) = 1. A more general analysis would lead to the 
result: 


3x + 3y + 3z — 4, 


(92) 


which is known as the Robinson damping theorem Q. The significance of this result is that while it 
is possible (for example, by changing the field gradient in the dipoles) to ‘shift’ the radiation damping 
between the different degrees of freedom, the overall amount of damping is fixed. In a planar storage 
ring, for example, one can reduce the horizontal damping time, but only at the expense of increasing the 
longitudinal damping time. 


In a typical storage ring, the dispersion in the dipoles is small compared to the bending radius of 
the dipoles, that is: 

n„ 

(93) 


^« 1 . 


P 


Then, if there is no quadrupole component in the dipoles (so that fci = 0 in the dipoles), comparing (35) 

(94) 


and (|62|) leads to: 


h 
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in which case: 


jx - 1, (95) 

jz - 2. (96) 

The horizontal damping time is approximately equal to the vertical damping time; the longitudinal damp¬ 
ing time is about half the vertical damping time. Typical values for the damping times in medium energy 
synchrotron light sources are some tens of milliseconds, or a few thousand turns. 


2.5 Quantum excitation 

So far, we have assumed a purely classical model for the radiation, in which energy can be radiated 
in arbitrarily small amounts. From the expressions for the evolutions of the emittances (301, (681 and 
(901, we see that if radiation was a purely classical process, the emittances would damp towards zero. 
However, quantum effects mean that radiation is emitted in discrete units (photons). As we shall see, 
this induces some ‘noise’ on the beam, known as quantum excitation, the effect of which is to increase 
the emittance. The beam in an electron (or positron) storage ring will eventually reach an equilibrium 
distribution determined by a balance between the radiation damping and the quantum excitation. In the 
remainder of this section, we shall derive expressions for the rate of quantum excitation and for the 
equilibrium emittances in an electron storage ring. 


In deriving the equation of motion (591 for the action of a particle emitting synchrotron radiation, 
we made the (classical) approximation that in a time interval dt, the momentum dP of the radiation 
emitted goes to zero as dt goes to zero. In reality, emission of radiation is quantized, so we are prevented 
from taking the limit dP —0. The equation of motion for the action (561 should then be written: 


~ I N{u)udu + 


Wl 

PqC 


W2 

w 


f 


N{u) du, 


(97) 


where N{u) is the number of photons emitted per unit time in the energy range from u to u + du. 
The first term on the right hand side of (97 1 just gives the same radiation damping as in the classical 
approximation; the second term is an excitation term that we previously neglected. 


To find an explicif expression for fhe rafe of change of fhe acfion in ferms of fhe beam and laffice 
paramefers, we need fo find expressions for fhe infegrals f N(u) u du and f N{u) u^ du. The required 
expressions can be found from fhe specfral disfribufion of synchrofron radiafion from a dipole magnef, 
which is anofher resulf fhaf we quofe from synchrofron radiafion fheory. The spectral distribution of 
radiation from a dipole magnet is given by Q : 


d^ 

d'd 



(98) 


where dV / dJ) is the energy radiated per unit time per unit frequency range, and = w/tCc is the radiation 
frequency uj divided by the critical frequency Uc'. 


LUc = 


3 7^c 

2T’ 


(99) 


P.y is the total energy radiated per unit time ( [sT] ), and is a modified Bessel function. Since fhe 

energy of a photon of frequency uj isu = huj, if follows fhaf: 


Using (981 and (1001, we find: 


1 dV 

N(u) du = - -— dJ}. 

rujj du 



N{u) udu = Pj, 


( 100 ) 

( 101 ) 
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and: 


( 102 ) 


Cq is a constant given by: 


noo -p 

/ N{u)u^du = 2Cq-i‘^—P^. 
lo P 


^ 55 h 

32\/3?7ic' 

-1-13, 


(103) 


For electrons (or positrons) Cq ^ 3.832 x 10 m. 


The next step is to substitute for the integrals in (971 from (1011 and (102 1 , substitute for wi and 


W 2 from (54 1 and (551, and average over the circumference of the ring. This gives an expression for the 


evolution of the horizontal action (for x <C and px <C Ppx)' 

d^x 2 2 2 -^5 

dt Tx JxTx I2 

where the fifth synchrotron radiation integral Is is given by: 


The T-L function {'Hx) is given by: 




Tdx — IxVx ‘^cy.xPxPpx T f^xP'. 


,2 

’px ■ 


(104) 


(105) 


(106) 


The damping time and horizontal damping partition number are given, as before, by ( [65| ) and 
Note that the excitation term is independent of the emittance: the quantum excitation does not simply 
modify the damping time, but leads to a non-zero equilibrium emittance. The equilibrium emittance eo 
is determined by the condition: 

d£r 

= 0, (107) 

£x=S0 


dt 


and is given by: 


£0 = CqY 


I 5 

jxl2 


(108) 


Note that eo is determined by the beam energy, the lattice functions (Courant-Snyder parameters and 
dispersion) in the dipoles, and the bending radius in the dipoles. We shall discuss how the design of the 
lattice affects the value of (and hence, the equilibrium horizontal emittance) in Section The equi¬ 
librium horizontal emittance ( 108| ) determined by radiation is sometimes called the natural emittance of 
the lattice, since it includes only the most fundamental effects that contribute to the emittance: radiation 
damping and quantum excitation. Other phenomena (such as impedance or scattering effects) can lead to 
some increase in the equilibrium emittance actually achieved in a storage ring, compared to the natural 
emittance. Typically, third generation synchrotron light sources have natural emittances of order of a few 
nanometres. With beta functions of a few metres, this implies horizontal beam sizes of tens of microns 
(in the absence of dispersion). 

In many storage rings, the vertical dispersion in the absence of alignment, steering and coupling 
errors is zero, so that Py = 0. However, the equilibrium vertical emittance is larger than zero, because 
the vertical opening angle of the radiation excites some vertical betatron oscillations. The fundamental 
lower limit on the vertical emittance, from the opening angle of the synchrotron radiation, is given by |[9||: 


13 Cn 


Sti — 


fj^ 

^^jyhJ \P^\ 


ds. 


(109) 


In most storage rings, this is an extremely small value, typically four orders of magnitude smaller than 
the natural (horizontal) emittance. In practice, the vertical emittance is dominated by magnet alignment 


17 
















Fig. 6: Change in longitudinal phase space variables for a particle emitting a photon carrying energy u. As a result 
of the photon emission, there is a change in amplitude of the synchrotron oscillations (represented by the ellipses) 
performed by the particle as it moves around the storage ring. 


errors. Storage rings typically operate with a vertical emittance that is of order 1% of the horizontal 
emittance, but many can achieve emittance ratios somewhat smaller than this. We shall discuss the 
vertical emittance in more detail in Section |4l 

Quantum effects excite longitudinal emittance as well as transverse emittance. Consider a particle 
with longitudinal co-ordinate z and energy deviation 5, which emits a photon of energy u (see Fig. [^. 
The co-ordinate and energy deviation after emission of the photon are given by: 


6' = Sq sin 0' = So sin 0 — —, 


z = 


Therefore: 


Averaging over the bunch gives: 


where: 


x> a' ^P^ 
—^On cos 0 = 


OJs 


OJs 


Eo 
do cos 0. 


— '^0 ~ 2(^0 sin 0 -|- 


Act? = 


Eo 

y) 


-^0 


2El- 


1 


= (5^) = 2 (^ 0 )- 


Including radiation damping, the energy spread evolves as: 


( 110 ) 

( 111 ) 

( 112 ) 

(113) 

(114) 


da'j 1 1 
dt 2 Eq Co 


dC 



duN{u)u^ - a'j, 

Ty 


(115) 


where we h ave averaged the radiation effects around the ring by integrating over the circumference. 
Using (102 1 for J iV {u)u‘^ du, we find: 


dt jzTzh 

The equilibrium energy spread is given by da'j/dt = 0: 



(116) 


^lo = 


Jl. 

jzh’ 


(117) 


18 














where the third synchrotron radiation integral I 3 is defined: 


h 



( 118 ) 


The equilibrium energy spread aso determined by radiation effects is often referred to as the natural 
energy spread, since collective effects can often lead to an increase in the energy spread with increasing 
bunch charge. Note that the natural energy spread is determined essentially by the beam energy and by 
the bending radii of the dipoles; rather counterintuitively, it does not depend on the rf parameters (either 
the voltage or the frequency). On the other hand, the bunch length does have a dependence on the rf. 
The ratio of the bunch length to the energy spread in a matched distribution (i.e. a distribution that 
is unchanged after one complete revolution around the ring) can be determined from the shape of the 
ellipse in longitudinal phase space followed by a particle obeying the longitudinal equations of motion 
([72]) and (f/^. Neglecting radiation effects (which can be assumed to be small) the result is: 


OpC 

(Tz = -^as- 

Us 


(119) 


We can increase the synchrotron frequency Ug, and hence reduce the bunch length, by increasing the rf 
voltage, or by increasing the rf frequency. 


2.6 Summary of radiation damping and quantum excitation 

To summarise, including the effects of radiation damping and quantum excitation, the emittances (in 
each of the three degrees of freedom) evolve with time as: 


s{t) 


e{t = 0 ) exp 



+ £(f 


00) 



( 120 ) 


where e{t = 0 ) is the initial emittance (for example, of a beam as it is injected into the storage ring), 
and e{t = 00 ) is the equilibrium emittance determined by the balance between radiation damping and 
quantum excitation. The damping times are given by: 


3oi 

II 

II 

N 

II 

bO 

( 121 ) 

where the damping partition numbers 

are given by: 


jx — i- ~ 

h . . h 

1^’ 7y — 1 ) ~ 1^' 

( 122 ) 

The energy loss per turn is given by: 


(123) 

where for electrons (or positrons) C-^ : 

8.846 X 10“^ m/GeV^. The natural emittance is: 



eo = Cqy^^, 

Jxh 

(124) 

where for electrons (or positrons) Cq 

Ri 3.832 X 10 m. The natural rms energy spread and bunch 

length are given by: 

2 r> ^,2 -^3 

Jzh 

(125) 


CXpC 

(Tz = - (TS- 

Us 

(126) 
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The momentum compaction factor is: 


(127) 


_ h 

Co' 

The synchrotron frequency and synchronous phase are given by: 

qVrl Wrf 
Eq To 

qVrf 

Finally, the synchrotron radiation integrals are: 



h 

h 

h 

h 

h 



P 

I 

J \P\' 


— { ^ + 2ki\ ds 


ds, 


ki = 


dB„ 


Pq dx 


'Hx = 


— ‘2‘O'xVxVpx + (^xVpx- 


(128) 

(129) 

(130) 

(131) 

(132) 

(133) 

(134) 


3 Equilibrium emittance and storage ring lattice design 


In this section, we shall derive expressions for the natural emittance in four types of lattices: FODO, 
double bend achromat (DBA), multi-bend achromat (including the triple bend achromat) and theoretical 
minimum emittance (TME) lattices. We shall also consider how the emittance of an achromat may be 
reduced by ‘detuning’ the lattice from the strict achromat conditions. 


Recall that the natural emittance in a storage ring is given by (108 1 : 


£0 = 


h 

jxh' 


(135) 


where Cq is a physical constant, 7 is the relativistic factor, jx is the horizontal damping partition number, 
and Is and I 2 are synchrotron radiation integrals. Note that jx, h and I 2 are all fixed by the layout of the 
lattice and the optics, and are independent of the beam energy. In most storage rings, if the bends have no 
quadrupole component, the damping partition number jx ^ 1- In that case, to find the natural emittance 
we just need to evaluate the two synchrotron radiation integrals I 2 and If we know the strength and 
length of all the dipoles in the lattice, it is straightforward to calculate I 2 . For example, if all the bends 
are identical, then in a complete ring (total bending angle = 27r): 


h 



B ds 2 'kB cB 

_= _ ~ Ott_ 

{Bp) p {Bp) E/q' 


(136) 


where E is the beam energy, and q is the particle charge. Evaluating Is is more complicated: it depends 
on the lattice functions. 


3.1 FODO lattice 

Eet us consider the case of a EODO lattice. The lattice functions in a typical EODO cell are shown 
in Eig. To simplify the system, we use the following approximations. Eirst, we assume that the 
quadrupoles can be represented by thin lenses. Second, we assume that the space between the quadrupoles 
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Fig. 7: Lattice functions in a FODO cell. Top; Courant-Snyder parameters and dispersion. Bottom: % function. 
In this case, the phase advance is 90°, the dipoles are 1.5 m long and have bending angle 27r/32. Notice that the 
value of the T-L function is constant except in the dipoles: this is a general property of this function. 


is completely filled by the dipoles. This is clearly not a realistic assumption, but it does allow us to de¬ 
rive some useful (and reasonably accurate) formulae. With these approximations, the lattice functions 
(Courant-Snyder parameters and dispersion) are completely determined by the focal length / of the 
quadrupoles and the bending radius p and length L of the dipoles, and can be calculated using standard 
techniques. 

Suppose that i^ceii is the transfer matrix for the horizontal motion in one complete periodic cell of 
a lattice. i?ceii may be constructed by multiplying the transfer matrices R for individual components in 
the beam line. For example, for a thin quadrupole of focal length /: 



0 

1 
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For a dipole of bending radius p and length L, the transfer matrix is: 


( cos 4 p sin 4 

l-^L £ 

__sm- cos- 


( 138 ) 


The Courant-Snyder parameters at any point in the beam line can be found first by multiplying the 
transfer matrices R for the individual components to give the transfer matrix i?ceii for the periodic cell 
starting from the chosen point, and then writing the complete transfer matrix in the form: 


( cos Px + Oix sin Px (3x sin px 

■^cell — I • 

y —7a; Sin ^ 2 ^ cos /ix — sm //a: 

where px is the phase advance. The dispersion describes the periodic trajectory of an (off-energy) particle 
through a periodic cell, and can be found at any point by solving the condition: 


(139) 


Vx 

Ppx 


- 7?^ 

~ “cell 


Px 

Ppx 


+ dl 


cell ’ 


(140) 


where is a matrix representing the first order terms in the map (for a complete cell) for the dispersion, 
and d£^^ is a vector representing the zeroth order terms. The map for a complete cell is found, as usual, 
by composing the maps for individual elements. For a quadrupole, the map for the dispersion is the same 
as the map for the dynamical variables; for a dipole, there are additional zeroth order terms: 


Px 

Vpx 


= R, 


dip 


so+^ 


Vx 

Vpx 



p{l — cos 

sin — 
p 



(141) 


Using the above results, we find that in terms of /, p and L, the horizontal beta function at the 
horizontally focusing quadrupole in a FODO cell is given by: 


Afp sin 9{2f cos 6* + p sin 9) 

716/^ — [p^ — (4/2 -|- p2 j cQg 2P]2 ’ 


(142) 


where 9 = L/p is the bending angle of a single dipole. The dispersion at a horizontally focusing 
quadrupole is given by: 


2 /p(2/ + ptan^) 

4/2 + p2 


(143) 


By symmetry, at the centre of a quadrupole, ax = Ppx = 0. Given the lattice functions at any point 
in the lattice, we can evolve the functions through the lattice, using the transfer matrices R. For the 
Courant-Snyder parameters: 


A(si) = RA{so)R^, 


(144) 


where R = R{si; sq) is the transfer matrix from sq to si, R^ is the transpose of R, and: 


A = 


Px 

Cix 



(145) 


The dispersion can be evolved (over a distance L, with constant bending radius p) using (141 1 . 


We now have all the information we need to find an expression for Is in the FODO cell. However, 
the algebra is rather formidable. The result is most easily expressed as a power series in the dipole 
bending angle, 0: 




JX-I 


- {j-A + OiO* 


(146) 
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Fig. 8: Ratio of synchrotron radiation integrals I^jli in a FODO cell, as a function of the phase advance. The 
black line shows the exact value, while the red line shows the value calculated using the approximation (1481. 


For small 0, the expression for I ^/12 can be written: 



This can be further simplified if p S> 2/ (which is often the case): 

h A \ 

up) p3’ 

and still further simplified if 4f ^ L (which is less oflen fhe case): 


(147) 


(148) 


h ' 


(149) 


The rafio /12 is ploffed for a FODO cell as a function of fhe phase advance in Fig. Making fhe 
approximafion jx ~ 1 (since we assume fhaf fhere is no quadrupole componenf in fhe dipole), and 
writing p = L/6, we have: 




(150) 


Nofice how fhe emiffance scales wifh fhe beam and laffice parameters. The emiffance is propor¬ 
tional fo fhe square of fhe energy and fo fhe cube of fhe bending angle. Increasing fhe number of cells 
in a complete circular laffice reduces fhe bending angle of each dipole, and reduces fhe emiffance. The 
emiffance is proporfional fo fhe cube of fhe quadrupole focal lengfh: sfronger focusing resulfs in lower 
emiffance. Finally, fhe emiffance is inversely proportional fo fhe cube of fhe cell lengfh. 

The phase advance in a FODO cell is given by: 


P 

cos Px — f 2 ■ 


(151) 


This means fhaf a sfable lattice musf have: 


/ 1 
- > -. 
L - 2 
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( 152 ) 









In the limiting case, Hx 
gives: 


180°, and / has the minimum value / = L/2. Using the approximation (1501 


and so the minimum emittance in a FODO lattice is expected to be: 


^0,FODO,min 




(153) 


However, as we increase the focusing strength, the approximations we used to obtain the simple expres¬ 
sion for £0 start to break down. From the exact formula for / 5 /I 2 as a function of the phase advance, we 
find (by numerical means) that there is a minimum in the natural emittance at pLx ~ 137° « 0.38 x 27r rad 
(see Fig.[^. The minimum value of the natural emittance in a FODO lattice is given by: 

eo,FODO,min ~ 1.2Cq7^0^. (154) 


As an example, consider a storage ring with 16 FODO cells ( 32 dip oles), 90° phase advance per 
cell (/ = and with a stored beam energy of 2GeV. Using (1501 we estimate that such a ring 

would have a natural emittance of around 125 nm. Many modem applications (including synchrotron 
light sources) demand emittances smaller than this by one or two orders of magnitude. This raises 
the question of how we might design a lattice with a smaller natural emittance. Looking at the lattice 
functions in a FODO lattice (Fig.[^ provides a clue. The dispersion function, which is directly related to 
the effect of quantum excitation on the horizontal emittance, is non-zero throughout the cell. If we can 
design a lattice where the dispersion vanishes at the entrance of a dipole, then we might hope to reduce 
the average value of the 'K function in the dipoles, thereby reducing Is and the value of the natural 
emittance. It is indeed possible to design a cell with two dipoles, in which the dispersion vanishes at the 
entrance of the first dipole and at the exit of the second dipole: such a cell is known as a Chasman-Green 
cell flOl, or a double bend achromat (DBA). 


3.2 Double bend achromat lattice 

To calculate the natural emittance in a DBA lattice, let us begin by considering the conditions for zero 
dispersion at the start and the exit of a unit cell. Assume that the dispersion is zero at the start of the 
cell. We place a quadrupole midway between the dipoles, to reverse the gradient of the dispersion. By 
symmetry, the dispersion at the exit of the cell will then also be zero. In the thin lens approximation, the 
required strength of the quadrupole between the dipoles can be determined from: 


V~i// 1/V Vpx ) V ~ ^ / V ~Vpx ) 


(155) 


Hence the central quadrupole must have focal length: 


Vx 

‘^Pjpx 


(156) 


The actual value of the dispersion (and its gradient) is determined by the dipole bending angle 6, the 
bending radius p, and the drift length Ldrift^ 


rjx = /9(1 - COS0)-h Ldriftsin 6 », (157) 

Ppx = sin 6 ». (158) 

To complete the DBA cell, we need to include some additional quadrupoles in the zero-dispersion 
region to control the horizontal and vertical beta functions. To correct the chromaticity, sextupoles are 
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Fig. 9: Lattice functions in a DBA cell. Top; Courant-Snyder parameters and dispersion. Bottom: % function. 
The horizontal beta and alpha functions at the entrance of the first dipole have values fix = 2.08 m and ax = 2.47. 
These are different from the ‘ideal’ values for low emittance in this case, of fix — 2.33 m and Ux = 3.87. The 
lattice functions are detuned from their ideal values in order to satisfy a range of constraints (such as maximum 
values of the beta functions, magnet strengths and chromaticity). The detuning results in this case in an increase in 
the natural emittance by a factor of 1.8. 


included between the dipoles, where the dispersion is non-zero. The lattice functions in an example 
DBA cell are shown in Fig.[^ To get some idea whether this style of lattice likely to have a lower natural 
emittance than a FODO lattice, we can inspect the T-L function. Comparing Figs. and we see that 
the T-L function is much smaller in the DBA lattice than in the FODO lattice. Note that we use the same 
dipoles (bending angle and length) in both cases. 

Let us calculate the minimum natural emittance of a DBA lattice, for given bending radius p and 
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bending angle 0 in the dipoles. To do this, we need to calculate the minimum value of: 

j-L nj 

h= -^ds (159) 

Jo 

in one dipole (of length L), subject to the constraints: 


Px,0 — Ppx,0 — 0 ) 


(160) 


where rjxfl and rjpx^ are the dispersion and gradient of the dispersion at the entrance of the dipole. 
We know how the dispersion and the Courant-Snyder parameters evolve through the dipole, so we can 
calculate for one dipole, for given initial values of the Courant-Snyder parameters Q;a:,o and /3xfl- 
Then, we have to minimise the value of Is with respect to axfi and Pxfl- Again, the algebra is rather 
formidable, and the full expression for is not especially enlightening: therefore, we just quote the 
signihcant results. We hnd that, for given p and 9 and with the constraints ( |160 l the minimum value of 
Is is given by: 

= ^- + 0(06). ( 161 ) 

4V15 P 

This minimum occurs for values of the Courant-Snyder parameters at the entrance to the dipole given 
by: 


[dx,o — 


L + o{e^), 


— \/l5 + O(0^), 

where L = is the length of a dipole. Since we know that I 2 in a single dipole is given by: 


1 


9 


d2= -^ds= 

Jo P P 

we can now write down an expression for the minimum emittance in a DBA lattice: 


^OjDBA.min — Cq'J 


2 -^5,min 

jxl2 


1 


4v/l5 




(162) 

(163) 

(164) 

(165) 


The approximation is valid for small 9. Note that we have again assumed that, since there is no 
quadrupole component in the dipole, jx ^ 

Compare the expression ( 165| ) for the minimum emittance in a DBA lattice, with the expression 
( 154| ) for the minimum emittance in a FODO lattice. We see that in both cases (FODO and DBA), the 
emittance scales with the square of the beam energy, and with the cube of the bending angle. However, 
the emittance in a DBA lattice is smaller than that in a FODO lattice (for given energy and dipole bending 
angle) by a factor 4:y/l5 « 15.5. 

This is a signihcant improvement; however, there is still the possibility of reducing the natural 
emittance (for a given beam energy and number of cells) even further. For a DBA lattice, we imposed 
constraints ( 160| ) on the dispersion at the entrance of the hrst dipole in a lattice cell. To reach a lower 
emittance, we can consider relaxing these constraints. 

3.3 Theoretical minimum emittance lattice 

To derive the conditions for a theoretical minimum emittance (TME) lattice, we write down an expression 
for: 

Ux 


h — 


ds, 


(166) 
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with arbitrary dispersion 'qxfl, r]px,o and Courant-Snyder parameters axfi and /3xfi in a dipole with given 
bending radius p and angle 6 (and length L = p9). Then, we minimise Is with respect to r]x,o, flpxfl, Oix,o 
and Pxfl- The result is (IT): 

£0,TME,mm ~ 12-\/l5 (16V) 

The minimum emittance is obtained with dispersion at the entrance to the dipole given by: 


Vx,o = ]:L9 + 0{e^), (168) 

b 

Vpx,o = -^ + 0(9^}, (169) 

and with Courant-Snyder functions at the entrance: 

/3xfi = -^L + 0{9^), (170) 

a:,,o = Vi5 + 0{9^). (171) 


The dispersion and beta function reach minimum values in the centre of the dipole: 


Vx,mm = ^ ^ 


(172) 

(173) 


By symmetry, we can consider a single TME cell to contain a single dipole, rather than a pair of dipoles 
as was necessary for the DBA cell. Outside the dipole, the dispersion is relatively large. This is not ideal 
for a light source, since insertion devices at locations with large dispersion will blow up the emittance. 
If insertion devices are required, then it is possible to break the symmetry of the lattice to include zero- 
dispersion straights: for example, the ring could have a race-track footprint, with arcs constructed from 
TME cells. 


Examples of the lattice functions (and Ti function) in a TME cell are shown in Eig. 10 Note that 
the Ti function in the dipole in the TME cell is significantly lower than for EODO or DBA cells using 
similar dipoles (Eigs.j^anc 


3.4 Practical constraints on lattice optics 

The results we have derived for the natural emittance in EODO, DBA and TME lattices have been for 
‘ideal’ lattices that perfectly achieve the stated conditions in each case. In practice, lattices rarely, if ever, 
achieve the ideal conditions. In particular, the beta function in an achromat is usually not optimal for 
low emittance; and it is difficult to tune the dispersion for the ideal TME conditions. The main reasons 
for this are: first, beam dynamics issues (relating, for example, to nonlinear dynamics and collective 
effects) often impose a variety of strong constraints on the design; and second, optimising the lattice 
functions while respecting all the various constraints can require complex configurations of quadrupoles. 
A particularly challenging constraint on design of a low-emittance lattice is the dynamic aperture. Stor¬ 
age rings require a large dynamic aperture in order to achieve good injection efficiency and good beam 
lifefime. However, low emittance lattices generally need low dispersion and beta functions, and hence 
require strong quadrupoles. As a result, the chromaticity can be large, and must be corrected using strong 
sextupoles. Strong sextupoles lead to highly nonlinear motion and a limited dynamic aperture: the tra¬ 
jectories of particles at even quite modest betatron amplitudes or energy deviations can become unstable, 
resulting in short beam lifetime. 

Eattices composed of DBA cells have been a popular choice for third generation synchrotron light 
sources. The DBA structure provides a lower natural emittance than a EODO lattice with the same 
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Fig. 10: Lattice functions in a TME cell. Top; Courant-Snyder parameters and dispersion. Bottom; % function. 
The horizontal beta function and dispersion match the ‘ideal’ values for low emittance. 


number of dipoles, while the long, dispersion-free straight sections provide ideal locations for insertion 
devices such as undulators and wigglers. If an insertion device, such as an undulator or wiggler, is 
incorporated in a storage ring at a location with large dispersion, then the dipole fields in the device 
can make a significant contribution to the quantum excitation (Is). As a result, the insertion device can 
lead to an increase in the natural emittance of the storage ring. By using a DBA lattice, dispersion-free 
straights are naturally provided, in which undulators and wigglers can be located without blowing up 
the natural emittance. However, there is some tolerance. In many cases, it is possible to detune the 
lattice from the strict DBA conditions, thereby allowing some reduction in natural emittance at the cost 
of some dispersion in the straights. The insertion devices will then contribute to the quantum excitation; 
but depending on the lattice and the insertion devices, there may still be a net benefit. Some light sources 
that were originally designed with zero-dispersion straights take advantage of tuning flexibility to operate 
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with non-zero dispersion in the straights (see, for example, p^). This provides a lower natural emittance, 
and better output for users. 


3.5 Multi-bend achromats 

There are of course many options for the design of a storage ring lattice, beyond the FODO, DBA 
and TME cells we have discussed so far. For example, it is possible to combine the DBA and TME 
lattices, constructing an arc cell consisting of more than two dipoles. The dipoles at either end of the cell 
have zero dispersion (and gradient of the dispersion) at their outside faces, thus satisfying the achromat 
condition. Since the lattice functions are different in the central dipoles compared to the end dipoles, 
we have additional degrees of freedom we can use to minimise the quantum excitation. The result is a 
multi-bend achromat (MBA) that combines the benefits of a DBA lattice (with long straights providing 
good locations for insertion devices) and a TME lattice (providing the possibility of achieving lower 
emittance than in a DBA). 

In a MBA, it is possible to have cases where the end dipoles and central dipoles differ in the 
bend angle (i.e. length of dipole), and/or the bend radius (i.e. strength of dipole). For simplicity, let us 
consider the case where the dipoles all have the same bending radius (i.e. they all have the same field 
strength), but they vary in length. Assume that each arc cell has a fixed number M of dipoles, with 
average bending angle 9 = 27r/MAceiis- If the two outer dipoles have bending angle aO and the inner 
dipoles have bending angle bO, then the coefficients a and b satisfy: 

2a + (M - 2)b = M. (174) 


Eet us assume that the lattice functions (Courant-Snyder parameters and dispersion) in the outer dipoles 
are the same as in a DBA lattice, and in the inner dipoles are the same as in a TME lattice. Since the 
synchrotron radiation integrals are additive, for an M-bend achromat, we can write: 


2 (aO)^ {M - 2) {bO f _ 6a^ + (M - 2)b^ 9^ 
^ 4^/l5 P ^ 12\/l5 P ~ 12\/l5 P ’ 

nf) hB f) 

/ 2 ,ceii « 2—+ (M-2)-= (2a + (M-2)6)-. 

P P P 

Hence, in an M-bend achromat: 

4, cell ^ 1 f 6a^ + {M - 2)6"^ ^ ^3 

/2,ceii ~ 12\/l5 \2a + {M -2)b ) 


(175) 

(176) 


(177) 


Minimising the ratio I 5 /I 2 with respect to a gives: 


a 1 
b~ 

from which it follows that: 

+ {M -2)b^\ _ M + 1 
V 2a + {M-2)b ) ~ M -I 


(178) 


(179) 


The central bending magnets should be longer than the outer bending magnets by a factor Then, the 
minimum natural emittance in an M-bend achromat is given by: 


^0,MBA,min 


1 

12^15 


/M + l\ 
\M-l) 




(180) 


Note that 9 is the average bending angle per dipole. Although we derived (180) with the assumption of 
at least three dipoles (M > 2), the formula gives the correct result for a DBA in the case M = 2. Also, 
in the limit M —)■ 00 , we obtain the correct expression for the natural emittance in a TME lattice. 
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Table 1: Minimum natural emittance in different lattice styles for electron storage rings: for each lattice style, the 
minimum natural emittance is given by , where Cq « 3.832 x 10“^^ m, and 7 is the relativistic factor 

for the beam. The dipoles have length L and bending angle Q, and no quadrupole component. 


Lattice style 

T 

Conditions 

90° FODO 

2 V 2 

f = L/V2 

137° FODO 

1.2 

minimum emittance FODO 

DBA 

1 

4vT5 

Vx,o — — 0? /3x,o ^ o(x,o ^ V^15 

MBA 

1 

12x/l5 \M-lJ 

M dipoles (with same radius of curvature) per cell 

TME 

1 

12^/15 

r? ■ K. B • Ri ^ 

7a:,mm ~ 24 ~ 2 \/T 5 


Triple bend achromats have been used in light sources, including the ALS 1131 and the SLS 1141. 
Light sources based on cells with even larger numbers of bends per achromat are planned: see, for 
example, 0 As with double bend achromats, it is possible to obtain some reduction in the natural 
emittance of a triple (or higher) bend achromat by detuning the lattice from the strict achromat condition, 
allowing some dispersion to Teak’ into the straight sections. As long as the dispersion in the straights is 
not too large, there is a net benefit, despite some contribution to the emittance from quantum excitation 
in the insertion devices. 

As a final remark, we note that further flexibility to optimise the natural emittance can be provided 
by relaxing the constraint that the field strength in a dipole is constant along the length of the dipole. We 
expect an optimised design to have the strongest field at the centre of the dipole, where the dispersion 


can be minimised. For an example, see 1161. 


4 Vertical emittance generation, calculation and tuning 

In this section, we shall discuss how vertical emittance is generated by alignment and tuning errors, 
describe methods for calculating the vertical emittance in the presence of known errors, and discuss 
briefly how an operating storage ring can be tuned to minimise the vertical emittance (even when the 
alignment and tuning errors are not well known). 


Recall that the natural (horizontal) emittance in a storage ring is given by (108 1 : 

2 h 


£0 = C'gT" 


3x^2 


(181) 


If the horizontal and vertical motion are independent of each other (i.e. if there is no betatron coupling) 
then we can apply the same analysis to the vertical motion as we did to the horizontal. If we build a ring 
that is completely flat (i.e. no vertical bending), then there is no vertical dispersion, i.e. r\y = r/py = 0 at 
all locations around the ring. It follows that the vertical 'H. function 'Ky. 

~ “b ‘^^yVy'npy + Pylpy: ( 182 ) 

also vanishes around the entire ring, and that therefore the synchrotron radiation integral I^y will be zero. 
This implies that the vertical emittance will damp to zero. 

However, in deriving equation ( 181| ) for the natural emittance, we assumed that all photons were 
emitted directly along the instantaneous direction of motion of the electron. In fact, photons are emitted 
with a distribution having angular width I /7 about the direction of motion of the electron. This leads to 
some vertical ‘recoil’ that excites vertical betatron motion, resulting in a non-zero vertical emittance. A 
detailed analysis leads to the following formula for the fundamental lower limit on the vertical emittance 

&■ 


13 a 


-y,mm 


_ 

^^jyhJ IpP 


ds. 


(183) 
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To estimate a typical value for the lower limit on the vertical emittance, let us write equation (1831 in the 
approximate form: 


-y,min 


^jyh 


J\P? ^ 


(184) 


where {j3y) is the average vertical beta function around the ring. Using some typical values {{fiy) = 20 m, 
jz = 2, jy = I, as = 10“^, 7 = 6000), we find: 


^y,min ~ 0.3 pm. 


(185) 


The lowest vertical emittance achieved so far in a storage ring is around a picometer, several times 
larger than the fundamental lower limit (see, for example, 181). In practice, vertical emittance in a 
(nominally planar) storage ring is dominated by two effects: residual vertical dispersion, which couples 
longitudinal and vertical motion; and betatron coupling, which couples horizontal and vertical motion. 
The dominant causes of residual vertical dispersion and betatron coupling are magnet alignment errors, 
in particular: tilts of the dipoles around the beam axis; vertical alignment errors on the quadrupoles; tilts 
of the quadrupoles around the beam axis; and vertical alignment errors of the sextupoles. Let us consider 
these errors in a little more detail. 


Steering errors lead to a distortion of the closed orbit, which generates vertical dispersion and 
(through vertical offsets of the beam in the sextupoles) betatron coupling. A vertical steering error may 
be generated by rotation of a dipole, so that the field is not exactly vertical, or by vertical misalignment 
of a quadrupole, so that there is a horizontal magnetic field af fhe location of the reference trajectory. 

Coupling errors lead to a transfer of horizontal betatron motion and dispersion into the vertical 
plane: in both cases, the result is an increase in vertical emittance. Coupling may result from rotation 
of a quadrupole, so that the field contains a skew component. When particles pass through a skew 
quadrupole, they receive a vertical kick that depends on their horizontal offset. As a result, quantum 
excitation of the horizontal emittance feeds into the vertical plane. 

A vertical beam offset in a sextupole has the same effect as a skew quadrupole. To understand 
this, recall that a sextupole field is given by: 


Bz, = {Bp)k2xy, (186) 

By = ]j^{Bp)k2[x^-y^) ■ (187) 

A vertical offset can be represented by the transformation y ^ y + Ay: 


Bx 1 -^ {Bp)k2xy + {Bp)k2Ayx, (188) 

By ^ ^{Bp)k2 (x^ - y^) - {Bp)k2Ayy - ^k2Ay‘^. (189) 


The terms in (188 1 and (1891 that are first order in Ay constitute a skew quadrupole of strength {Bp)k 2 Ay. 


When designing and building a storage ring, we need to know how accurately the magnets must be 
aligned, to keep the vertical emittance below some specified limit. Although beam-based tuning methods 
also normally have to be applied, the ultimate emittance achieved after machine tuning does depend on 
the accuracy with which the initial alignment is performed. It is therefore useful to have expressions 
that relate the closed orbit distortion, vertical dispersion, betatron coupling and (ultimately) the vertical 
emittance, to the alignment errors on the magnets. 


4.1 Closed orbit distortion 

Let us begin by considering the closed orbit distortion. In terms of the action-angle variables, we can 
write the coordinate and momentum of a particle at any point: 

y = ^/2PyJy cos (t)y, (190) 
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Fig. 11: Closed orbit distortion from a thin dipole kick in a synchrotron storage ring. If the co-ordinate and 
momentum of a particle on the closed orbit immediately after the dipole kick are {yo,Pyo), then after nearly one 
complete turn, just before the dipole kick the co-ordinate and momentum of the particle are {yo,Pyo ~ A0). The 
dipole kick then puts the particle back onto the closed orbit. 


Py 



(sin + ay cos (py) . 


(191) 


Suppose there is a steering error at some location s = sq which leads to an instantaneous change (i.e. a 
‘kick’) A9 in the vertical momentum. After one complete turn of the storage ring, starting from imme¬ 
diately after sq, the trajectory of a particle will close on itself if: 


^‘^PyoJyO COS (pyl — '^/‘^PyoJyO COS <pyO, 


I 2JyO 
PyO 


(sin ayo cos (py\) — 


I 2JyO 

f^yO 


sin (pyo ayo cos (pyo) AO . 


(192) 


(193) 


where cpyi = (pyQ + 2Tnyy, and Vy = //y/27r is the vertical tune (see Fig. 111. Solving equations (192 1 
and (|193|) for the action and angle at sq: 


JyO 

4^y0 


(3yo^e^ 

n • 9 ’ 

S Sin TTUy 

TTUy. 


(194) 

(195) 


Note that if the tune is an integer, there is no solution for the closed orbit: even the smallest steering error 
will kick the beam out of the ring. From (1951, we can write the coordinate for the closed orbit at any 
point in the ring: 


( \ y/ Pyi^o)Py{s) ^ ^ 

Ucois) = ^ . _ A0cOS['KVy + pLy[s\ Sq)) , 


2 sin TTVr, 


(196) 


y 


where /iy(s; sq) is the phase advance from sq to s. 

In general, there will be many steering errors distributed around a storage ring. The closed orbit 
can be found by summing the effects of all the steering errors: 


Vcois) 


ypyjs) 

2 sin TTKy 



+ p.y{s\ s')) els'. 


(197) 


It is often helpful to be able to estimate the size of the closed orbit distortion that may be expected 
from random quadrupole misalignments of a given magnitude. We can derive an expression for this from 
equation ( 197| ). For a quadrupole of integrated focusing strength kiL, vertically misaligned from the 
reference trajectory by AY, the steering is: 


AO = {kiL)AY. 


(198) 
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Estimated sensitivity: 19.1651 (simulation), 15.23 (analytical) 



RMS quadrupole vertical misalignment [jim] 


Fig. 12: Simulation of closed orbit distortion resulting from quadrupole alignment errors in a storage ring GD- 
Each circle shows the mean of the rms orbit distortion from 100 different sets (seeds) of random alignment errors on 
the quadrupoles; the error bars show the range covered by 90% of the seeds. The solid red line shows a linear fit to 
the circles; the broken red line shows an analytical estimate of the orbit distortion based on the known quadrupole 
strengths and lattice functions, using equation (|199|l. 


Squaring equation (197 1 , then averaging over many seeds of random alignment errors, we find: 

'yl{s)\_ (Ay2) 


I3y{s) 


) Sin TTU. 


E 

quads 


PyihLY 


(199) 


In performing the average, we assume that the alignments of different quadrupoles are not correlated in 
any way. 

The ratio between the closed orbit rms and the magnet misalignment rms is sometimes known as 
the orbit amplification factor. Values for the orbit amplification factor are typically in the range from 
10 to about 100. Of course, the amplification factor is a statistical quantity: the actual rms of the orbit 
distortion depends on the particular set of alignment errors present. 

In the context of low-emittance storage rings, vertical closed orbit errors are of concern for two 
reasons. First, vertical steering generates vertical dispersion, which is a source of vertical emittance. 
Second, vertical orbit errors contribute to vertical beam offset in the sextupoles, which effectively gen¬ 
erates skew quadrupole fields, which in turn lead to betatron coupling. We have seen how to analyse the 
beam dynamics to understand the closed orbit distortion that arises from quadrupole alignment errors of 
a given magnitude. Our goal is to relate quantities such as orbit distortion, vertical dispersion, coupling, 
and vertical emittance, to the alignment errors on the magnets. We continue with betatron coupling. 


4.2 Betatron coupling 

Betatron coupling describes the effects that can arise when the vertical motion of a particle depends on 
its horizontal motion, and vice-versa. Betatron coupling can arise (for example) from skew quadrupoles 
and solenoids. 

In a storage ring, skew quadrupole fields ofen arise from quadrupole fills, and from vertical align¬ 
ment errors on sextupoles. A full treatment of betatron coupling can become quite complex, and there 
are many different formalisms that can be used. However, it is possible to use a simplified model to 
derive approximate expressions the equilibrium emittances in the presence of coupling. The procedure 
is as follows. First, we write down the equations of motion for a single particle in a beamline containing 
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coupling. Then, we look for a ‘steady state’ solution to the equations of motion, in which the horizontal 
and vertical actions are each constants of the motion. Finally, we assume that the actions in the steady 
state solution correspond to the equilibrium emittances (since e = (J)), and that the sum of the horizon¬ 
tal and vertical emittances is equal to the natural emittance of the ‘ideal’ lattice (i.e. the natural emittance 
of the lattice in the absence of errors). This procedure can give some useful results, but because of the 
approximations involved, the formulae are not always very accurate. 

We will use Hamiltonian mechanics. In this formalism, the equations of motion for the action- 
angle variables (with path length s as the independent variable) are derived from the Hamiltonian: 


H — H e/x, (jiy, Jy] s), 

using Hamilton’s equations: 


dJx 

dH 

ds 

d4>x 

dJy 

dH 

ds 

d(j)y 

dff^x 

dH 

ds 

dJx' 

d(f>y 

dH 

ds 

dJy 


For a particle moving along a linear, uncoupled beamline, the Hamiltonian is: 

TT _ Jy 

Px Py 


( 200 ) 

( 201 ) 

( 202 ) 

(203) 

(204) 

(205) 


The first step is to derive an appropriate form for the Hamiltonian in a storage ring with skew 
quadrupole perturbations. In Cartesian variables, the equations of motion in a skew quadrupole can be 
written: 


where: 


dpx 

ds 

= ksV, 

(206) 

dpy 

ds 

= ksx, 

(207) 

dx 

d,s 

= Px, 

(208) 

dv 

ds 

= Py^ 

(209) 

ks = 

1 dBx 

Bp dx 

(210) 


These equations can be derived from the Hamiltonian: 


H = -pi + -Py - ksxy. 


( 211 ) 


We are interested in the case where there are skew quadrupoles distributed around a storage ring. 
The ‘focusing’ effect of a skew quadrupole is represented by a term in the Hamiltonian: 

kgxy = 2ksy/I3xl3yy^JxJy cos cos (j)y. (212) 
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This implies that the Hamiltonian for a beam line with distributed skew quadrupoles can be written: 


H = ^ ^ - 2ks{s)^/^^^/JJ^cos4>xCos(|)y. (213) 

Px Py 

The beta functions and the skew quadrupole strength are functions of the position s. This makes it 
difficult to solve the equations of motion exactly. Therefore, we simplify the problem by ‘averaging’ the 
Hamiltonian: 

H = OJxJx + ^yjy — JxJy COS 4>x COS (l)y. (214) 

Here, ujx, ^y are the phase advances per unit length of the beam line, given by: 


1 /■'^o ds 

C'o Jo Px,y 


(215) 


where Co is the circumference of the ring, k is a constant that characterises the coupling strength. For 
reasons that will become clear shortly, we re-write the coupling term, to put the Hamiltonian in the form: 


H — UJxdx T ^y^y \l Jxdy COs((/)3; 4^y^ y/ Jxdy COs((/)3; 4“ (j^y^• 


(216) 


The constants k± represent the skew quadrupole strength averaged around the ring. However, we need 
to take into account that the kick from a skew quadrupole depends on the betatron phase. Thus, we write: 


= 


1 

a 





(217) 


where fix and are the betatron phase advances from the start of the ring. 

Now suppose that /t_ » /t+. (This can occur, for example, if lOx ^ u)y, in which case all the 
contributions to R- from the skew quadrupole perturbations will add together in phase.) Then, we can 
simplify things further by dropping the term in /t+ from the Hamiltonian: 

H — uJxJx 4“ ^y^y y/ JxJy cos((^^ ^y')' (718) 


We can now write down the equations of motion: 


dJx 

ds 

dH 

d4>x 

— K— Jx 

Jy sm{(j)x - 

4’y)i 

dJy 

ds 

dH 

d(t>y 

— JxJy 

sm{4>x 

~ 4>y ): 


dH 

K- 


COs(^(px 


ds 

dJx 

u;. + —^ 

/ 

— (f)y) 

^ 1 

11 

dH 

dJy 

R- 

jh. 
1 Jx 

COsi^Cpx 

— 4>y) 


(219) 

( 220 ) 

( 221 ) 

( 222 ) 


Even after all the simplifications we have made, the equations of motion are still rather difficult to 
solve. Fortunately, however, we do not require the general solution. In fact, we are only interested in the 
properties of some special cases. First of all, we note that from (|219|) and (|220[): 


ddx dJu „ 
- 7 ^ + -?^ =0 
ds ds 


(223) 


and therefore the sum of the actions Jx + Jy is constant. Going further, we notice that if 4>x = <t>y’ then 
the rate of change of each action falls to zero. This implies that if we can find a solution to the equations 
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Fig. 13: Variation of the ‘fixed point’ actions (|226|) and \221) as a function of the strength of the coupling reso- 


of motion with (px 
and: 


then: 


(py for all s, then the actions will remain constant. In fact, we find that if (px = (py^ 


(i(px d(py 

ds ds ’ 


(224) 


jy _ y^1 + «-/Aw2 - 1 

+ R'ilAoj‘^ + 1 


(225) 


where Atu = ojx — u)y. If we further use Jx + Jy = Jo, where Jq is a constant, then we have a solution 
to the equations of motion in which the actions are constant, and given by: 


1 


1 


Jx — 9 1 H- - - 

^ + K?./Aa;2^ 


Jv = ^ I 1 - 


1 


+ K^/Aa;2 


Jo, 


Jo- 


(226) 


(227) 


Note the behaviour, shown in Fig. [T^ of the fixed actions as we vary fhe ‘coupling sfrengfh’ k_ 
and fhe befafron funes (befafron frequencies). The fixed acfions are well-separafed for k_ <C Acj, buf 
bofh approach fhe value Jo/2 for Auj. The condifion af which fhe funes are equal (or differ by an 

exacf infeger) is known as fhe difference coupling resonance. 

Recall fhaf fhe emiffance may be defined as fhe befafron acfion averaged over all particles in fhe 

beam: 

Bx = {Jx), and £y = (Jy). (228) 

Now, synchrofron radiafion will damp fhe beam fowards an equilibrium disfribufion. In fhis equilibrium, 
we expecf fhe befafron acfions of fhe particles fo change only slowly, i.e. on fhe fimescale of fhe radiafion 
damping, which is much longer fhan fhe fimescale of fhe befafron motion. In thaf case, fhe acfions of 
mosf parficles musf be in fhe correcf rafio for a fixed-poinf solufion fo fhe equations of mofion. Then, if 
we assume fhat £x + £y = eo> where Eq is fhe nafural emiffance of fhe sforage ring, we musf have for fhe 
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Fig, 14: Effect of a single skew quadrupole (at a location with zero dispersion) on the vertical emittance 
synchrotron storage ring, as a function of the difference in the betatron tunes. The circles show the results 
computation using Chao’s method |211; the red line shows an analytical estimate using equation (230l. 


in a 
of a 


equilibrium emittances: 


£x — 1 + 


£y — 


= 1 - 


Y^l + K^jAuj"^ 
1 

+ k?_/Au;2 


2 ’ 


2 ■ 


(229) 


(230) 


As an illustration, we can plot the vertical emittance as a function of the ‘tune split’ At/, in a model 
of the ILC damping rings, with a single skew quadrupole (located at a point of zero dispersion, so as not 


to couple horizontal dispersion into the vertical plane). The result is shown in Fig. 14 The tunes are 
controlled by adjusting the regular (normal) quadrupoles in the lattice. The simulation results are based 
on emittance calculation using Chao’s method, which we shall discuss later. 

The presence of skew quadrupole errors in a storage ring affects the betatron tunes. To estimate 


the size of the effect, we use the Hamiltonian (2181. If we consider a particle close to the fixed point 
solution, we can assume that cpx = so that the Hamiltonian becomes: 


H — ^xJx T \/ JxJy' 


(231) 


The normal modes describe motion that is periodic with a single well-defined frequency. In the absence 
of coupling, the transverse normal modes correspond to motion in just the horizontal or vertical plane. 
When coupling is present, the normal modes involve a combination of horizontal and vertical motion. 


Let us write the Hamiltonian (231 1 in the form: 


H = 


( y/Jx ^/Jy ) A ^ ^ , 


where: 


A = 


-iff- 


-ift- 

LOti 


(232) 


(233) 


The normal modes can be constructed from the eigenvectors of the matrix A, and the frequency of each 
mode is given by the corresponding eigenvalue. From the eigenvalues of A, we find thaf the normal mode 
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Fig. 15: Effect of a single skew quadrupole (at a location with zero dispersion) on the measured betatron tunes in a 
synchrotron storage ring, as a function of the difference in the betatron tunes in the absence of the skew quadrupole. 
The circles show the results of a computation of the eigenvalues of the single-turn transfer matrix; the solid lines 
show an analytical estimate using equation ( |235| l. The minimum difference between the measured tunes gives the 
coupling strength in the ring. 


frequencies are: 



Hence, the tunes i/± are given (in terms of the tunes Ux and i/y in the absence of errors) by: 


(234) 



+ r'j/ 




(235) 


where, from (217 1 , k 


(Co/27r)K_ is given by: 


(236) 

The dependence of the tunes on the coupling strength provides a useful method for measuring the cou¬ 
pling strength k in a real lattice. The procedure is simple: a quadrupole (or combination of quadrupoles) 
is used to change the tunes, and then the tunes are recorded and plotted as a function of quadrupole 
strength. The minimum separation between the measured tunes gives the coupling strength. An exam¬ 
ple (from simulation) is shown in Fig. [T^ Of course, this procedure does not identify the source of 
the coupling, or provide very much information as to an optimal correction (beyond the strength of a 
skew quadrupole that may be required to achieve the correction, assuming that the skew quadrupole is 
at the correct phase in the lattice). However, the technique can be useful to characterise the effect of a 
correction that may need to be applied in several iterations. 


= 


1 

2vr Jo 


Major sources of coupling in storage rings are quadrupole tilts and sextupole alignment. Using 
the theory just outlined, we can estimate the alignment tolerances on these magnets, for given optics and 
specified vertical emittance. Starting with equation (236), we first take the modulus squared, and then 
use (for a sextupole with vertical alignment error Ay 5 ) kg = k 2 ^Ys and (for a quadrupole with tilt error 
A0 q) kg = A:iA0q. Assuming that there are no correlations between the errors, we find: 




My{k2l? + 

sexts 


(Ae^) 

dvr 2 


I3xl3y{k,lf, 

quads 


(237) 


where (k^) represenfs fhe mean value of fhe square of fhe coupling sfrengfh over a large number of sefs 
of random errors. Nofe fhaf AI 5 is fhe beam offsel from fhe cenfre of a sexfupole: fhis includes fhe 
effecfs of closed orbif disforfion. 
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4.3 Vertical dispersion 

Vertical emittance is generated by vertical dispersion as well as by betatron coupling. Vertical dispersion 
is in turn generated by vertical closed orbit distortion (vertical steering), and coupling of horizontal 
dispersion into the vertical plane by skew quadrupole fields. Our goal now is to estimate the amount of 
vertical dispersion generated from magnet alignment errors; we can then estimate the contribution to the 
vertical emittance. 

The equation of motion for the vertical co-ordinate for a particle with momentum P is: 


d^y ^ = In 

ds^ (Bp) P ■ 


(238) 


For small energy deviation 5, P is related to the reference momentum Pq by: 


(l + <5)Po. 

We can write for the horizontal field (to first order in the derivatives): 

„ „ dB^ dB^ 

Bx ~ Bqx + y —p. -h X ——. 

oy ox 

If we consider a particle following an off-momentum closed orbit, so that: 


(239) 


(240) 


y = Vyd, (241) 

X = rjxd, (242) 


then, combining the above equations, we find to first order in 6: 

d'^Py 


ds^ 


kiPy ~ k^g -|- k\sPx- 


(243) 


Equation (2431 gives the ‘equation of motion’ for the dispersion. It is similar to the equation of 
motion for the closed orbit: 

- kiyco ~ -kos + kigXco- (244) 

We can therefore immediately generalise the relationship ( |199| ) between the closed orbit and the quadrupole 
misalignments, to find for the dispersion: 


pI 

/dy 




Q! 


) sm TTZ/., 




quads 


> sm TTZ/., 


^ plPy{kiLf + '^plPy{k2Lf 


quads 


> sm 


y sexts 


(245) 

Here, we assume that the skew dipole terms /cqs come from vertical alignment errors on the quadrupoles 
with mean square (Alg), and that the skew quadrupoles ku come from tilts on the quadrupoles with 
mean square (A0 q) and from vertical alignment errors on the sextupoles, with mean square (AVj). We 
assume that all alignment errors are uncorrelated. 

The final step is to relate the vertical dispersion to the vertical emittance. This is not too difficult. 
First, we can apply the formula ( |108| ) for the natural (horizontal) emittance to the vertical emittance: 

(246) 




jyh 


where jy is the vertical damping partition number (usually, jy = 1), and the synchrotron radiation 
integrals are given by: 

r P.. 

(247) 


- y |^|3 
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and: 


(248) 


h = 


ds. 


The vertical V. function is: 


~ 'yy^y ‘^OiyVyVpy + PyVpy 


(249) 


If the vertical dispersion is generated randomly, then we can assume that it will not be correlated 
with the curvature l/p of the reference trajectory. (This is not the case for the horizontal dispersion!) 
Then, we can write: 

I^y ~ {Tdy) ^ ds = 


(250) 


Hence, for the vertical emittance: 


Ey ^ Cy-i^l(H.y)- 


h 


jyh 


It is convenient to use (1171 for the mean square energy spread, to give: 

Ey K, ^{'Hy)crg. 

Jy 

Now, note the similarity between the action: 

2Jy = + 2ayypy + /SyP^, 

and the Ti function: 


This implies that we can write: 
and hence: 


^y ~ ^y^y ‘^ay'TjyPpy + PyPpy. 


Tjy - jdyldy cos (j^yy^ 




= 


(251) 

(252) 

(253) 

(254) 

(255) 

(256) 


Combining equations (252i and (2561 gives a useful (approximate) relationship, between the vertical 
dispersion and the vertical emittance: 


dy \/^y / 


(257) 


Equation (245 I tells us how the vertical dispersion depends on the magnet alignment, and equation 
( 257| ) tells us how the vertical emittance depends on the vertical dispersion. Simply combining these two 
equations gives us an equation for the contribution of the vertical dispersion to the emittance, in terms of 
the magnet alignment errors. 

It should be remembered that the total vertical emittance is found by adding together the contri¬ 
butions from betatron coupling (equations ( 230| ) and (237 1 ) and vertical dispersion (equations (2451 and 
(257])). All these expressions involve significant approximations. However, they can give results that 
agree reasonably well with more reliable methods: an example is shown in Fig.[T^ 
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Fig. 16: Simulation of vertical emittance resulting from sextupole alignment errors in a storage ring m Each 
circle shows the mean of the normalised vertical emittance {'ySy) from 100 different sets (seeds) of random align¬ 
ment errors on the sextupoles; the error bars show the range covered by 90% of the seeds. The solid red line shows 
a quadratic fit to the circles; the broken red line (very close to the solid red line) shows an analytical estimate of the 
emittance based on the known sextupole strengths and lattice functions, using equations ( |230| l, ( |237| l to estimate 
the coupling contribution, and (|245[) and (257|i to estimate the dispersion contribution. 


4.4 Accurate computation of emittance 

The formulae we have derived so far are useful for developing a ‘feel’ for how the vertical emittance 
depends on magnet alignment errors, and for making rough estimates of the sensitivity to particular types 
of error. For detailed studies, including modelling and simulations, we need more accurate formulae for 
computing the vertical emittance in a storage ring with a given set of alignment errors. The calculations 
involved then become quite complex, and need to be solved using a computer. 

There are three methods commonly used for computing the equilibrium emittances in complex 
lattices with known errors. First, there is a technique based on the usual formulae for the emittances 
expressed in terms of the radiation integrals, but generalised to the normal modes (see, for example, 
g). Second, there is Chao’s method | [^ , which involves integrating the eigenvectors of the single¬ 
turn transfer matrix around the circumference of the ring. Finally, there is the ‘envelope’ method p^ , in 
which the second order moments of the equilibrium beam distribution are first computed from the single- 
turn transfer map (including radiation damping and quantum excitation); then the emittances are obtained 
from the matrix describing the beam distribution. We shall discuss briefly each of these techniques in 
turn. 


First, we consider the method for computing the equilibrium emittances based on normal mode 
analysis. Let us assume that we have a lattice code that will compute the symplectic single-turn trans¬ 
fer matrix at any point in a given lattice. In general, the transfer matrix will have non-zero terms off 
the block-diagonals: these terms represent coupling between the horizontal, vertical, and longitudinal 
motion. The expression (1811 we derived for the natural emittance assumed no betatron coupling, and 
that the coupling between the horizontal and longitudinal motion was relatively weak. However, we can 
generalise the formula to the case that betatron coupling is present. We still need to assume that the 
longitudinal motion is weakly coupled to each of the transverse degrees of freedom (i.e. the horizontal 
and vertical motion). In that case, we can consider separately the 4x4 single-turn transfer matrix 
describing the transverse motion: 

R±_ • 

• R\\ 


R = 


(258) 
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i?ll is a 2 X 2 matrix describing the longitudinal motion, and we assume we can neglect the terms repre¬ 
sented by the bullets (•). 

Now we look for a transformation, represented by a 4 x 4 matrix V, that puts R± into block- 
diagonal form, i.e. that ‘decouples’ the transverse motion: 


and i?n are 2 x 2 matrices describing betatron motion in a coordinate system in which the motion 
appears uncoupled. There are various recipes for constructing the decoupling transformation V (which 
is not unique): see, for example ||2^24|. Having obtained the matrices describing the uncoupled motion, 
we can derive the Courant-Snyder parameters for the normal mode motion in the usual way. For example, 
we can write: 

cos /ill + cqi sin /in Ai sin /in 

—711 sin /in cos /in — an sin /in 


Rn = 


(260) 


and similarly for mode I. We can also obtain the normal mode dispersion functions, by applying the 
transformation H to a vector constructed from the dispersion functions in the original Cartesian co¬ 
ordinates. Then, we can construct the R function for each mode; for example: 


Rn = Tii^n + 2an?7ii7p,ii + Aif/pji- 


(261) 


Finally, we can write for the mode II emittance: 


eii = Cq'y'^ 


I5,11 

h — h,n ’ 


(262) 


and similarly for mode I. 

For many storage rings, equation (2621 works well, and gives an accurate result. However, if there 
is strong coupling between the transverse and the longitudinal motion (which can happen, for example, 
for large values of the synchrotron tune), then the approximations needed to derive equation ( 262) start 
to break down. 


As an alternative to the normal mode analysis, we can consider Chao’s method | |2T| for computing 
the emittances, which provides a formula that can be expressed in a convenient form, though it is not 
always easy to apply. It is again based on the single-turn transfer matrix, but it is more accurate than the 
‘decoupling’ method, since it uses the full 6 x 6 transfer matrix, and does not assume weak coupling 
between the longitudinal and transverse motion. We do not explain here the physics behind the formula, 
but simply quote the result: 


ek = CL^ 
cak 


\Ek,{s)\^ 

|p(s)|3 


ds, 


(263) 


where A: = I, II, III is an index that specifies a particular degree of freedom, the eigenvalues of the 
single-turn matrix including radiation damping are is the fifth component of the A:* 

eigenvector of the symplectic single-turn matrix, and: 


55 Teh 
48^3 rn ’ 


(264) 


where Tc is the classical radius and m the mass of the particles in the beam. 


Finally, we mention the envelope method |22|. Like Chao’s method, it gives accurate results for 
the emittances even if there is strong coupling between all three degrees of freedom. The envelope 
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method is based on finding the equilibrium beam distribution described by the Sigma matrix: 
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(X2) 

{xpx) 
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{PxV) 

{PxPy) 

(PxZ) 
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(yx) 
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(y^) 


{PyX) 

(PyPx) 

{Pyy) 

(pD 

{Pyz) 
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{zpx) 

{zy) 

{zPy) 

{z^) 

{z5) 

V 

{6x) 

i^Px) 

i^y) 

i^Py) 

{6z) 

{5^) 


This is a symmetric matrix, constructed from the second order moments of all possible combinations 
of the dynamical variables. For simplicity, we assume in what follows that the first order moments 
are all zero, i.e. that the closed orbit lies along the reference trajectory. However, the method is easily 
generalised to include cases where there is closed orbit distortion. In the absence of coupling, the Sigma 
matrix will be block diagonal. We are interested in the more general case, where coupling is present. 

Under a single turn around an accelerator, S transforms as: 

S I—)> RTiR^ + D, (266) 


where R is the single-turn transfer matrix (including radiation damping) and U is a constant matrix rep¬ 
resenting the effects of quantum excitation. From knowledge of the properties of synchrotron radiation, 
we can compute the matrices R and D for a given lattice design: this will be discussed further below, 
where we shall give explicit expressions for the transfer matrices in a dipole, including radiation effects. 

The equilibrium distribution Sgq has the property: 

-f D. (267) 


For given R and D, we can solve equation (267 1 to find Sgq, and then from Sgq we can find the invariant 
emittances, i.e. the conserved quantities under symplectic transport. For any beam distribution S, the 
invariant emittances are given by: 


eigenvalues(SS') = ±isk, 


(268) 


where S is the antisymmetric block-diagonal matrix Q. To see that this is the case, consider the (simpler) 
case of motion in one degree of freedom. The Sigma matrix in this case is: 


/ (x^) (xpa;) ^ ^ Px -Ox 

V (pI) ) V 


In one degree of freedom, the matrix corresponding to Q is: 


S = 


0 

-1 


1 

0 


(269) 


(270) 


Then, the eigenvalues of US' are ztie^. Now, we can show that (under certain assumptions) the emittance 
is conserved as a bunch is transported along a beam line. In any number of degrees of freedom, the linear 
transformation in phase space co-ordinates of a particle in the bunch between two points in the beam line 
can be represented by a matrix R: 

X I—)■ Rx, (271) 

where x is a vector whose components are the phase space variables x,. 

Now consider how the Sigma matrix transforms. The Sigma matrix can be written as the product 
of the phase-space co-ordinates averaged over the bunch: 


Sjj — (xjXj), 


(272) 
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where Y^ij is the {i,j) component of the Sigma matrix, and the Xi are the dynamical variables. The 
brackets (•) indicate an average over all particles in the bunch. Then, using (2711, it follows that under a 
transformation R of the dynamical variables, the Sigma matrix transforms as: 


S KYR^ 


(273) 


Since S' is a constant matrix, it immediately follows that: 

YS ^ RYR^S. (274) 

Then, using the fact that R is symplectic ([^, we have: 

YS^ RYSR-^. (275) 


This is a similarity transformation of 115: the eigenvalues of any matrix are conserved under a similarity 
transformation. Therefore, since the eigenvalues of S5 give the emittance of the bunch, it follows that 
the emittances are conserved under linear, symplectic transport. 

This argument applies for any number of degrees of freedom. We define the matrix 5 in three 
degrees of freedom by (|^. The six eigenvalues of YS are then ±fefc, where k is an index ranging over 
the different degrees of freedom. The quantities are all conserved under linear, symplectic transport. 
Even if, as is generally the case, the Sigma matrix is not block-diagonal (i.e. if there is coupling present), 
then we can still find three invariant emittances using this method, without any modification. 

Neglecting radiaton, if 7? is a (symplectic) matrix that represents the linear single-turn transforma¬ 
tion at some point in a storage ring, then an invariant or ‘matched’ distribution is one that satisfies: 

S RYR^ = Y. (276) 


In general, all fhe particles in the bunch change position in phase space after one turn around the ring: 
but for a matched distribution, the second order moments remain the same. Although this condition 
determines the lattice functions (which can be found from the eigenvectors of YS), it is not sufficient to 
determine the emittances. In other words, the matched distribution condition determines the shape of the 
bunch, but not the size of the bunch. This makes sense: after all, in a proton storage ring, we can have 
a matched bunch of any emittance. However, in an electron storage ring, we know that radiation effects 
will damp the emittances to some equilibrium values. We shall now show how to apply the concept 
of a matched distribution, when radiation effects are included, to find the equilibrium emittances in an 
electron storage ring. 

To account for radiation effects in an electron storage ring, we must make two modifications to the 
single-turn transformation. First, the matrix R will no longer be symplectic: this accounts for radiation 
damping. Second, as well as first order terms in the transformation (represented by the matrix R), there 
will be zeroth order terms: these will correspond to the quantum excitation. The condition for a matched 
distribution should then be written: 

S = RYR'^ + D, (277) 


where R and D are constant, non-symplectic matrices that represent the first order and zeroth order 
terms in the single-turn transformation, respectively. Equation (277 1 is sufficient to determine the Sigma 
matrix uniquely - in other words, using just this equation (with known R and D) we can find the bunch 
emittances and the matched lattice functions. 


The envelope method for finding the equilibrium emittances in a storage ring then consists of 
three steps. First, we need to find the first order terms R and zeroth order terms D in the single-turn 
transformation: 


S RYR'^ + D. 


(278) 
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In the second step, we use the matching condition {111 i to determine the Sigma matrix. Then, in the 
third and final step, we find the equilibrium emittances from the eigenvalues of T,S. 


Strictly speaking, since R is not symplectic, the emittances are not invariant as the bunch moves 
around the ring. Therefore, we may expect to find a different emittance at each point around the ring. 
However, if radiation effects are fairly small, then the variations in the emittances will also be small. 


The transfer matrices R and D for an entire ring can be constructed from the transfer matrices for 
individual components in the ring. As an example, we shall consider a thin ‘slice’ of a dipole. This is an 
important case, since in most storage rings, radiation effects are significant only in dipoles. Furthermore, 
complete dipoles can be constructed by composing the maps for a number of slices. Hence, once we 
have a map for a thin slice of dipole, and knowing the usual (symplectic) transfer maps for drift spaces, 
quadrupoles and rf cavities, we will be able to construct the map for one complete turn of a storage ring, 
starting at any point. 


Recall that the transformations for the phase space variables in the emission of radiation carrying 
momentum dP are: 
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I-)- 

X, 


(279) 



/ dP^ 
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(281) 
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(282) 
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V-p. 

]Py^ 
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z, 


(283) 



. dP 



s 


Pq 


(284) 


where Pq is the reference momentum. In general, dP is a function of the co-ordinates. To find the trans¬ 
formation matrices R and D, we find an explicit expression for dP/P q, and then write down the above 
transformations to first order. For an ultra-relativistic particle, the momentum lost through radiation can 
be expressed in terms of the synchrotron radiation power Py (energy loss per unit time): 


d^ 


Pi , 
-^dt 


Eq 



ds 

c 


(285) 


where p is the radius of curvature of the reference trajectory. The radiation power Py is given by ( [3T] ). In 
general, the dipole may have a quadrupole gradient, so the field is: 


B = Bq + Bix. (286) 

Also, the particle may have some energy deviation, so the total energy is: 

E = Eo{l + 6). (287) 

Substituting these expressions, we find (after some manipulation): 

Pi = (;^ + 2^1 (1 + (288) 

where ki is the normalised quadrupole gradient in the dipole: 

A:i = J-Pi. (289) 

Pq 

Hence, the normalised momentum loss may be written: 


45 



Expanding to first order in the phase space variables, we can write: 
dP C^Ej^ , / 1 


^^0 


C-y El , 
zvr 


2tt \p 


+ 2k\ - 


PJ P 


^xds + 2^^6ds + 0(x2) + 0(6^). 
2tt p^ 


( 291 ) 


Given the expression (2911 for dP/P q, the transformations (279l-(284i become (to first order in 
the dynamical variables): 
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The first order terms give the components of Rdip{ds), the transfer matrix for a thin slice (length ds) of 
a dipole. There is a zeroth order term in the map for the dynamical variables that will contribute to (the 
( 6 , 6 ) component of) Ddip{ds), which contains the zeroth order terms in the transformation of the Sigma 
matrix through a thin slice of a dipole. Since the ( 6 , 6 ) component of Ddipj ds) represents the quantity 
(A(5^), the contribution to this component from the zeroth order term in (297 1 will be second order in ds. 
We still have to take proper account of the quantum nature of the radiation. This will make an additional 
contribution to Ddip (ds). 

The zeroth order term in the map for the Sigma matrix is given by: 


[Z7dip(ds)]6g “ ^ ^ ^ 


(n 2 ) 


(298) 


where {u^) is the mean square of the photon energy. Using (102i, we find that, to zeroth order in the 
phase space variables: 

((S' 



277 p-^ 


(299) 


Note that this term is first order in ds, whereas the contribution to Ddip{ds) that we found previously 
was second order in ds. Hence, in the limit ds —)• 0, the latter contribution dominates over the previous 
contribution. 

Collecting the above results, and taking only dominant contributions in the limit ds —)• 0, we find 
that the radiation in a thin slice of dipole has an effect on the Sigma matrix that can be represented by: 


S I— )• Rdip{ds)TiRj:^p{ds) + Ddip{ds) 
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where: 
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and: 
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To construct the full single-tum transformation, we need to compose the maps for all the elements 
in the ring, including the radiation effects in the dipoles. It is straightforward to do this numerically using 
a computer. However, some care is needed in handling the D matrices. For example, given the Sigma 
matrix at a location sq, we find the Sigma matrix at a location si = sq + ds from: 


S(si) = R{si; so)T.{so)R^{si; So) + D{si;so). (303) 


Then the Sigma matrix at S 2 is given by: 

S(s2) = i?(s2;si)S(si)f?'^(s2;si)+ T)(s2;si) 

= R{s2-, so)T.{so)R^{s2; sq) + R{s2; si)D{si; so)R'^{s2; si) + D{s2; si). (304) 

Hence: 


R{s2-,So) 

= i?(s2;'Si)-R(si;so) 

(305) 

D{s2-,So) 

= R{s2; si)D{si] so)R'^is 2 ; si) + D{s 2 ] si). 

(306) 

Continuing the process, we find: 



R{Sn]So) 

= R(^Sn'i Sn—l)R{Sn—l'j Sn— 2 ) ' ' ' Risi] Sq) 

(307) 
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When composing the transfer maps for thin slices of a dipole, we have to remember to ‘interleave’ the 
radiation maps with the usual symplectic transport map for a thin slice of dipole. 


The next step in finding the equilibrium emittances is to solve the matching condition (277 1 to find 
fhe Sigma matrix for the equilibrium distribution. To do this (for given matrices R and D), we make use 
of the eigenvectors U of R, and the diagonal matrix A constructed from the eigenvalues of R: 


RU = AU. 


Defining S and D by: 


E = 

D = UDU^, 


fhe solution for the Sigma matrix can be written: 


Ejj — 
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(310) 

(311) 


(312) 


The above formulae enable us to find the matched (equilibrium) distribution S; the eigenvalues of SS" 
are then AiSk, where Sk are the emittances. 
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The envelope method makes explicit the fact that vertical emittance is generated by coupling 
between the vertical and longitudinal planes in regions where radiation is emitted (i.e. by vertical disper¬ 
sion in dipoles), and by coupling between the vertical and horizontal planes in regions where radiation 
is emitted (i.e. by betatron coupling in dipoles). Here, we need to be careful in the use of the term ‘cou¬ 
pling’. In this context, coupling means the presence of non-zero components off the block-diagonals 
in the single-turn matrix, R. Full characterisation of the coupling requires complete specification of all 
the components off the block-diagonals in R. Depending on these components, it is possible to have 
coupling in a storage ring, and not generate any vertical emittance. For example, one could construct 
a closed ‘coupling bump’ using sets of skew quadrupoles in a straight section in a storage ring. With 
proper care in the design, outside the region between the skew quadrupoles, the vertical motion can be 
completely decoupled from the horizontal and the longitudinal. Then, despite the presence of strong 
coupling in some parts of the storage ring, the equilibrium vertical emittance will come only from the 
opening angle of the cone describing the spatial distribution of the synchrotron radiation. 

Numerical computational procedures (such as the envelope method) for finding fhe equilibrium 
beam disfribufion in a sforage ring are imporfanf because fhey provide ways fo calculate fhe equilibrium 
emiffances in complex, coupled laffices. If is possible fo include ofher non-symplecfic effecls in fhe 
calculafion (such as, for example, infrabeam scaffering). 


4.5 Ultra-low emittance tuning 

Often, coupling comes from magnet alignment errors, which will not be known completely in an oper¬ 
ating machine. At the design stage, it is important to characterise the sensitivity of a lattice to magnet 
alignment errors, particularly regarding the vertical emittance. Being able to compute the beam emit- 
tances in a storage ring with coupling errors present makes it possible to simulate the effects of various 
types and sizes of alignment error - and then to optimise the lattice design to minimise the sensitivity to 
the likely errors. However, in practice, tuning a storage ring to achieve a vertical emittance of no more 
than a few picometres (which may be required for some applications) is a considerable challenge, even 
in a lattice designed so as to minimise the sensitivity to coupling errors. Accurate alignment (by survey) 
of the magnets is always the first step in achieving ultra-low emittances; but beam-based tuning methods 
will then also be needed. 

A variety of beam-based methods for tuning storage rings have been developed over the years. A 
typical procedure might look as follows: 


1 . 

2 . 

3. 


Align the magnets by a survey of the ring. Typically, quadrupoles need to be aligned to better than 
a few tens of microns, and sextupoles to better than a couple of hundred microns. 

Determine the positions of the BPMs relative to the quadrupoles. This is known as ‘beam-based 
alignment’ (BBA, see Fig. [TT] ), and can be achieved by steering the beam to a position in each 
quadrupole where changing the quadrupole strength has no effect on the orbit p5|. 


Correct the orbit (using steering magnets) so that it is as close as possible to the centres of the 
quadrupoles. 


4. Correct the vertical dispersion (using steering magnets and/or skew quadrupoles, and measuring 
the dispersion at the BPMs) as close to zero as possible. 

5. Correct the coupling, by adjusting skew quadrupoles so that an orbit ‘kick’ in one plane (from any 
orbit corrector) has no effect on the orbit in the other plane. 


Usually, the last three steps need to be iterated several (or even many) times. 

Results from the tuning procedure described above can be limited by systematic errors on the 
BPMs, which can affect dispersion and coupling measurements. A useful technique for overcoming 
such limitations is to apply Orbit Response Matrix (ORM) analysis |261. This can be used to determine 
a wide range of magnet and diagnostics parameters, including coupling errors and BPM tilts. Although 


48 





Fig. 17: Beam-based alignment in a quadmpole. If the beam passes off-axis through a quadrupole magnet, then 
varying the strength of the magnet changes the trajectory downstream of the magnet. A change in trajectory can be 
observed in a beam position monitor (BPM). One method of beam-based alignment consists of steering the beam 
(using upstream orbit corrector magnets) until changing the quadrupole strength has no effect on the beam position 
observed in the BPM. 


vertical emittances of order 1 pm have now been achieved (representing an emittance ratio of less than 
0 .1%), tuning an electron storage ring to operate in this regime still remains a challenging goal, requiring 
extensive work and application of a range of techniques to reduce errors. Even making measurements of 
emittances less than a few picometres is not straightfoward, and requires specialist instrumentation (see, 
for example, |[27|). 
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